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ON A KINETIC FITZHUGH-NAGUMO MODEL OF NEURONAL NETWORK 


S. MISCHLER, C. QUININAO, J. TOUBOUL 


Abstract. We investigate existence and uniqueness of solutions of a McKean-Vlasov evolution 
PDE representing the macroscopic behaviour of interacting Fitzhugh-Nagumo neurons. This 
equation is hypoelliptic, nonlocal and has unbounded coefficients. We prove existence of a solu¬ 
tion to the evolution equation and non trivial stationary solutions. Moreover, we demonstrate 
uniqueness of the stationary solution in the weakly nonlinear regime. Eventually, using a semi¬ 
group factorisation method, we show exponential nonlinear stability in the small connectivity 
regime. 
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1. Introduction 

This paper undertakes the analysis of the existence and uniqueness of solutions for a mean-field 
equation arising in the modeling of the macroscopic activity of the brain. This equation describes 
the large-scale dynamics of a model of the central nervous system, taking into account the fact 
that it is composed of a very large number of interconnected cells that manifest highly nonlinear 
dynamics and are subject to noise. Non-linearities in the intrinsic dynamics of individual cells are an 
essential element of the neural code. Indeed, nerve cells constantly regulate their electrical potential 
depending on the input they receive. This regulation results from intense ionic exchanges through 
the cellular membranes. The modeling of these dynamics led to the development of the celebrated 
Hodgkin-Huxley model m, a very precise description of ion exchanges through the membrane 
and their effects on the cell voltage. A simplification of this model conserving the most prominent 
aspects of the Hodgkin-Huxley model, the Fitzhugh-Nagumo (FhN) model [T?, 29], has gained the 
status of canonical model of excitable cells in neuroscience. This model constitutes a very good 
compromise between versatility and accuracy on the one hand, and relative mathematical simplicity 
on the other hand. It describes the evolution of the membrane potential v of the cell coupled to an 
auxiliary variable x , called the adaptation variable. Different neurons interact through synapses 
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that are either chemical or electrical. In the case of electrical synapses for instance, the evolution 
of the pair voltage-adaptation for a set of n neurons {{v\, x\), 1 < i < n} satisfy the equations: 


( 1 . 1 ) 


[ dv\ = (v\ (v\ - A) (1 - v l t ) - x\ + E" =1 - v 3 t ) + I t ) dt + adW l t 

I dx\ = \-ax\ + bv\) dt, 


where the cubic nonlinearity accounts for the cell excitability, I t is the input level, a and b are 
positive constants representing timescale and coupling between the two variables, and the processes 
{(Wl)t>Oi 1 < i < n} are independent Brownian motions accounting for the intrinsic noise at the 
level of each cell. In the sequel, for sake of simplicity, we assume that a 2 = 2 and It = Iq 6 R 
constant, but it is likely that some of our analysis can be extend to I t e L°°{ R + ) converging rapidly 
when t goes to infinity. The coefficients J,j represent the effect of the interconnection of cell j onto 
cell i. These coefficients are positive, and incorporate the information of the connectivity map. 
Under relatively weak assumptions on the distribution of these coefficients (see Appendix [A} , it 
is relatively classical to show that the system enjoys propagation of chaos property and finite sets 
of neurons converge in law towards a process whose density solves the McKean-Vlasov evolution 
PDE: 

(1.2) d t f = Q e [ c / f ]f:=d x (Af) + d v (B e (Sf)f) + % v f on(0,oo)xR 2 , 

(1.3) A = A(x,v) = ax - bv, B e (/ f ) = B(x, v ; e, <//), 

(1.4) B(x,v;s,j)=v(v-X)(v-l) + x-£(v-j) + I 0 , </> = ,/(/)= f vf(x,v)dvdx, 

J 

where e denotes the averaged value of the connectivity coefficients Jij and / = f(t,x,v ) > 0 is the 
density function of finding neurons with adaptation and voltage ( x,v ) € R 2 at time t > 0. The 
evolution equation (1.2) is complemented by an initial condition 

/(0,v) = /o(v)>0 in R 2 . 

Since the PDE can be written in divergence form, the initial normalization of the density is con¬ 
served. In particular, consistent with the derivation of the system, we have: 


/ f{t, x, v) dxdv = / fo(x, v) dxdv = 1, 
JR 2 


when /o is normalized. Moreover, the nonnegativity is also a classical result of this kind of equations 
(for a brief discussion see Section [3]), therefore we assume in the sequel that / is a probability 
density. 

From the mathematical viewpoint, this equation presents several interests. First, the system is 
not Hamiltonian and the dynamics may present several equilibria, therefore, methods involving a 
potential and its possible convexity may not be used. Second, intrinsic noise acts as a stochastic 
input only into the voltage variable (since it modifies the voltage through random fluctuations of 
the current), leaving the adaptation equation unchanged and yielding to a hypoelliptic equation. 
From the phenomenological viewpoint, this system is particularly rich. It shows a number of 
different regimes as parameters are varied, and in particular, as a function of the connectivity 
level: the system goes from a non-trivial stationary regime in which several stationary solutions 
may exist for strong coupling, to periodic solutions, and eventually to a unique stationary solution 
for weak coupling. This is illustrated in section [6j in particular, we present some numerical results 
of (1.11 for a large number of interacting neurons. 

In order to rigorously analyse equation (|1A 


we restrict ourself to the latter regime, and we shall 
demonstrate the existence, uniqueness and stability of solutions to the McKean-Vlasov equation 
in the limit of weak coupling. More precisely, we shall prove existence of solution and non trivial 


stationary solution to the evolution equation (1.2) without restriction on the connectivity coefficient 


e > 0, and next uniqueness of the stationary solution and its exponential NL stability in the small 
excitability regime. 
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1.1. Historical overview of macroscopic and kinetic models in neuroscience. As men¬ 
tioned above, the problem we study lies within a long tradition of works in the domain of the 
characterization of macroscopic behaviors in large neuronal networks. First efforts to describe the 
macroscopic activity of large neuron ensemble can be traced back to the work of Amari, Wilson 
and Cowan in the 1970s |2 BUSSES], where were introduced heuristically derived equations on the 
averaged membrane potential of a population of neurons. These models made the assumption that 
populations interact through a macroscopic variable, the averaged firing rate of the population, as¬ 
sumed to be a sigmoidal transform of the mean voltage. This model has been extremely successful 
in reproducing a number of macroscopic behaviors in the cortex, one of the most striking being 
related to pattern formation in the cortex associated to visual hallucinations m ( see also [?j f° r a 
recent review on the subject). The relatively simplicity and good agreement with neurological phe¬ 
nomena motivated to understand the relationship between the dynamics of individual cells activity 
and macroscopic models. This has been an important piece of work in the 1990s in the bio-physics 
community, using simplified (non-excitable) models and specific assumptions on the architecture 
of the network, including the assumption of sparse and balanced connectivity (the sum of all in¬ 
coming input vanishes). The sparse connectivity assumption was used by the authors to stated 
that the activity was uncorrelated and resulted in characterizing different neuronal states. 

Alternative approaches were also developed based on population density [T§J methods. These 
yield complex partial differential equations, that were reduced to a set of moment equations from 
which authors may deduce the behavior of the system. The validity of these moment reduction 
and their well-posedness is a complex issue debated in the literature, see e.g. |24| . A transition 
Markov two-states model governing the firing dynamics of the neurons in the network was recently 
introduced. In these models, the transition probability of the system, written through a master 
equation, is then handled using different physics techniques including van Karnpen expansions or 
path integral methods. This modeling recently gathered the interest of the community (see for 
example da E H3 EH)- 

The mathematical community also undertook the analysis of the problem since the beginning 
of this decade. In that domain, one can distinguish also two distinct approaches: on one side, the 
development of mathematical models for simplified or phenomenological neuronal models, and on 
the other side works on the precise neuronal models. The dynamics of solutions of macroscopic 
limits of phenomenological neuron models is much more developed. The characterization of the 
stationary (or periodic) solutions was done in a simplified model, the Wilson-Cowan system, which 
has the important advantage to yield a Gaussian solution whose mean and standard deviation sat¬ 
isfy a deterministic dynamical system that may be studied analytically j33( 1351 using the analysis 
of ordinary differential equations. Artificial spiking neuronal models representing the discontin¬ 
uous dynamics of the time to the next spike were analyzed in a number of situations, including 
construction of periodic solutions to the limit equation in the presence of delays [301 na EU in 
the same vein, an important result was demonstrated on integrate-and-fire models in the presence 
of noise and excitation: it was shown that too much excitation could prevent the existence of 
solutions for all times, as the firing rate blows up in finite time m- These approaches make use 
of functional analysis of PDEs and nonlocal age-structured type of equations. 


1.2. Organization of the paper. The paper is organized as follows. Section [5] summarizes 
our main results that are demonstrated in the rest of the paper. Section [3] is interested with 
the existence, uniqueness and a priori estimates on the solutions to the evolution equation, as 
well as, the existence of stationary solutions. The next sections prove the stability of the unique 
stationary solution. Our proof uses factorization of the linearized semigroup allowing to prove linear 
stability, which we complete in section [5] by an analysis of the nonlinear stability of the stationary 
solution. Along the way, a number of open problems were identified beyond the small connectivity 
regime treated here that we present in section [6] together with numerical simulations: we will 
observe that the stationary solution splits into two stable stationary solutions as connectivity is 
increased, and in an intermediate regime, periodic solutions emerge. Two appendices complete the 
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paper. Appendix |A| investigates the microscopic system and its convergence towards the mean-field 
equation (|1.2|) and Appendix [B| deals with the strict positivity of stationary solutions. 


2. Summary of the main results 

2.1. Functional spaces and norms. We start by introducing the functional framework in which 
we work throughout the paper. For any exponent p € [1, oo] and any nonnegative weight function 
oj, we denote by L p (w) the Lebesgue space L P (R 2 ; w dx dv) and for k 6 N the corresponding Sobolev 
spaces W k ’ p (R 2 -,codxdv). They are associated to the norms 

Wfhw = IMl>, I/II^.P( W ) = 1/IIL(.) + t \\ D x, V f\\ P Lp( u) - 

1=1 

For k > 1, we define the partial ^-derivative space W k (uj) by 

w v k ’ p (w) := {/ 6 W k -^^)-,Dlf e L p (co)}, 
and it is natural to associate them to the norm 

l/I^C) = l/lw-. P (.) + \\D k v f\\ P LHu y 

A particularly important space in our analysis, denoted by H 2 (u>), is 

H 2 (uj) = W 2,2 (uj) = {/ e H 1 (w) such that d 2 vv f e A 2 (w)}, 

together with the set of functions with finite entropy 

L 1 log L 1 := | / 6 L 1 ( K 2 ) such that / > 0 and < oo j, 

where we use the classical notation ■= f R2 f log/. Finally, for k > 0, let us introduce the 

exponential weight function: 

(2.1) m = e K ( M_1 ) w ith M := 1 + x 2 /2 + v 2 12. 

In the sequel, we will be brought to vary the constant k involved in the definition of m, therefore 
we introduce the shorthand rrii = e Ki ( M_1 ), AN. Unless otherwise specified, these sequences are 
constructed under the assumption that the sequence Ki is strictly increasing. 


2.2. Main results. We start by stating a result related to the well possedness of (1.2) and to the 
a priori bounds on the solution. Using classical theory of renormalized solutions, it is not hard to 
see that equation (1.21 has indeed weak solutions , which we naturally define as: 

Definition 2.1. Let /o be a normalized nonnegative function defined on R 2 such that ^(f 0 ) is 
well defined. We say that ft(x,v ) := (t,x,v) n- f(t,x,v ) is a weak solution to (1.2) if the following 
conditions are fulfilled: 

- f € G([0, oo); L 1 (M 2 )); 

- for almost any t> 0, / > 0 and 

/ f(t,x,v) dx dv = / /o(x, v) dxdv = 1; 

- for any p 6 G 1 ([0, oo); C™ (R 2 )) and any t > 0 it holds 

( 2 - 2 ) f pf t = f pfo + f f [d t p + d 2 v p-Ad x p-B e (f (f s ))d v p]f s . 

JK Z J 0 


Equipped with this definition we can state the 
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Theorem 2.2. For any /o e L 1 (M 2 ) nL 1 logL 1 nP(R 2 ), there exists a unique global weak solution 
ft to the FhN equation (1.2), that moreover satisfies 

(2-3) ll/tlU^M) < max(C' 0 , ||/o||l 1 (m))> 

and depends continuously in L 1 (M) to the initial datum. More precisely, if f n , o /o wi lfi(M) 
and J8£(f n ,o) ^ C then f n> t -*■ /t m L 1 (M) /or any later time t> 0. 

Furthermore, there exist two norms | ■ fl^i and || • ||^2 equivalent respectively to fl • ||#i( m ) and 
| • such that the following estimates hold true: 

(2-4) \\ft\\m(m) < max((7i, |/o||ii( m )), 

as well as 

(2-5) l/tl-Hi <max(C' 2 ,||/o||«i), 

and 

(2-6) WftWni - max(C 3 , ||/o||-hj)> 

where C\,C 2 ,C% are positive constants. 

The other two main results of the present work can be summarized in the following 

Theorem 2.3. For any £ > 0, there exists at least one stationary solution G e to the FhN statistical 
equation (1.2), that is 

G e € H 2 (m) nP(f 2 ), 0 = d x (AG e ) + d v {B s {/ Ge )G e ) + d 2 v G e in R 2 . 

R such that r}{e) -» 0 and such that any 

£->0 


(2.7) 

Moreover, there exists an increasing function rj :. 


solution to (2.7) satisfies 

II G - Go||i 2 (m) ^ 7 l( e )) 

where G o is the unique stationary solution corresponding to the case e = 0. 


Theorem 2.4. There exists e* > 0 such that, in the small connectivity regime e e (0,e*), the 
stationary solution is unique and exponentially stable. More precisely, there exist a* < 0 and 


r]*{e) :] 


, with rf (e)- *■ oo, such that if 


£-»-0 
rl 


/ 0 6 77 i (m)nP(R 2 ) and ||/o - G\ H i {m) < rf(e), 
then there exists C* = C*(fo,e*,e) > 0, such that 


\\ft-G\\ L 2 (m) <C* e° 


Vt > 0, 


where f t is the solution to (1.2) with initial condition /q. 


2.3. Other notations and definitions. We prepare to the demonstration of these results by 
introducing a few notations that will be used throughout the paper. For two given Banach spaces 
(e, HU) and (£, | ■ ||f ), we denote by £${E,£) the space of bounded linear operators from E to 
£ and we denote by fl ■ \\gg(E,£) the associated operator norm. The set of closed unbounded linear 
operators from E to £ with dense domain is denoted by c <ci’(E,£ ). In the special case when E = £, 
we simply write 38(E) = 38(E, E) and <tf(E) = ^(E, E). 

For a given aeK, we define the complex half plane 

A Q :={ze(C, Re(z) > a}. 

For a given Banach space X and A € 'ta(X) which generates a semigroup, we denote this associated 
semigroup by t > 0), by D( A) its domain, by N(A) its null space, by i?(A) its range, and 

by £(A) its spectrum. On the resolvent set p( A) = C\ £(A) we may define the resolvent operator 
K(A) by 


VzeC, 'R.a(z) (A - z) -1 . 
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Moreover, 6 J3{X) and has range equal to D( A). We recall that £ e E(A) is called an 

eigenvalue of A if iV (A - £) t {0}, and it called an isolated eigenvalue if there exists r > 0 such that 

E(A)n{zeC, |z-£| < r} = {£}. 

Since the notion of convolution of semigroups will be required, we recall it here. Let us consider 
some Banach spaces Ai, A 2 and X 3 and two given functions 

S 1 eL 1 ([0,o O y,&(X 1 ,X 2 )) and S 2 6 L\[0, oo);#(A 2 , A 3 )), 
one can define S 2 * S± 6 L 1 ([0, 00 ); 3§(X\, X 3 )) by 

(S 2 *Si)(t):= S 2 (t-s)S 1 (t)da, Vi> 0. 

Jo 

In the special case Si = S 2 and X± = A 2 = A 3 , is defined recursively by S^* 1 ^ = S and 

g(*n) _ g * g(*(n- 1 )) f or n> 1 Equipped with this definition, we state the 

Proposition 2.5. Let X. Y be two Banach spaces such that Y c A'. Let us consider S(t ) a 
continuous semigroup such that for all t > 0 

1 S(t)la ix) <C x e aU , X e {X,Y}, 

for some a* € K and positive constants Cx and Cy. If there exists © > 0 and Cx,y > 0 such that 
\\S(t)f\\ Y < Cx,Yt~ 6 e aU \\f\\x, V/eX,te(0,1], 
then, there exists neN, and a polynomial p n (t ) such that 
(2.8) \\S { * n) (t)f\\Y < Pn{t) e“’ 4 || f\\ x , V / e A, t > 0. 

In particular, for any a> a*, it holds 

l|S ( * n) (i)/b * C a , n e at \\f\\ x , V / e A, t > 0, 
for some positive constant C Q>n . 

This general result has been already established and used in [19j and [26], but we give an 
alternative, and somehow simpler, proof of it. 


Proof. Let us start by noticing that for X e {A, Y}, if 

(2-9) \\S ( * n) (t)f\\ x <pZ(t)e aU lfl x , Vt>0, 

for n e N and p*(£) a polynomial, then 


\S« n+ 1 »(t)f\\ x < [ t \\S(t-s)S(* n Hs)f\\ X ds<pZ + 1 (t)e a * t \\f\\x, 

Jo 

for pX+i = Cx Jq Pn ( s ) ds. So, by an immediate induction argument we get (2.9) for any n > 1 and 
Pn Y) ■ („_1)! • 

Let us now fix t e (0,1] and, without lost of generality, assume that 0 i N. In that case, if 
(2.10) \\S^ n \t)f\\ Y <C n t-^- n+ ^ e aU \\f\\x, Vie (0,1], 

for some n e N and C n a positive constant, then 


||S ( * (n+1)) (£)/||y 


< f t/2 \\S(t-s)S^ n \s)f\\ Y ds+ r \\S(t-s)S^ n \s)f\\ Y ds 

Jo Jt /2 

f' t//Q i t 

< / C x , Y {t- S )- @ e a *^\\S(* n \s)f\\xd S + / C Y e a ^ t - s) \\S^ n \s)f\\ Y ds 

Jo Jt /2 

r £/2 j,, n t * / \ 

< / C x , Y (t - s)~ e e a t Pn(s)\\f\\x ds + / C Y e a t C n s~ ( - e ~ n+1 ^\\f\\ x ds 

JO Jt/2 

C* C in rt/2 rt 

< , X,Y * e^lflx [ {t-s)- e s n - 1 ds + C Y C n e a * t \\f\\x f A (e -" +1) ds 

(n - 1)! Jo Jt /2 

< C n+ it< @ - n h aU \\f\\x, 
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for some C n+ 1 depending only on Cx,Cy,Cx,y and C n . Once again, by an induction argument, 
we get (pOOl). Moreover, as soon as 0 - n + 1 > 0, inequality p.8|) holds. 


Finally, to get the conclusion in the case t > 1, it suffices to notice that 

II'S'W/IIy ^ CyCyy(t-[tJ) 6 e“ ‘l/f; 


\\X i 


where [fj is the largest integer smaller than t. A similar argument that the one used for t € (0,1], 
allows us to find a polynomial p n such that (2.8) still holds when t > 1. □ 


Finally, we recall the abstract notion of hypodissipative operators: 

Definition 2.6. Considering a Banach space (X, || • ||a')> a real number a € R. and an operator 
A e 'if(X), (A -a) is said to be hypodissipative on X if there exists some norm |||•||| Y on X equivalent 
to the usual norm || ■ ||a such that 

V/eD(A), 3</>eF(/) such that Ae(</>, (A - a)/) < 0, 

where (•,•) is the duality bracket in X and X* and F(f ) c X* is the dual set of f defined by 

F(f ) = F Hx U) := {<t>* X*, (0,/) = \jf\fx = Mix.}. 

One classically sees (we refer to for example [19], Subsection 2.3]) that when A is the generator 
of a semigroup S\, for given a € ffi. and C > 0 constants, the following assertions are equivalent: 

(a) (A - a) is hypodissipative; 

(b) the semigroup satisfies the growth estimate ||S'a(0II^(.y) ^ Ce at , t> 0; 


3. Analysis of the nonlinear evolution equation 

This section is concerned with the analysis of the nonlinear evolution equation. We shall prove 
existence and uniqueness of solutions, and provide some a priori estimates on their behavior. 

Before going into further details, let us remark that for Jf fixed, the operator is linear 

and writes 

QeU] f = d x {Af) + d v {Bdj?) f) + dlj. 

In particular, for g € H%(m) we have 

fAQe[/]f)gdvdx = - ff (Ad x g + B e (f) d v g - d 2 vv g) dvdx, 
therefore, it is natural to define 

Q*e[/]g : = ~Ad x g - Bdf) d v g + dl v g. 


3.1. A priori bounds. We now fix £o > 0. The a priori estimates that follow are uniform in e in 
the bounded connectivity regime e e [0 ,£q)j he., they involve constants that do not depend on e. 


Lemma 3.1. For f t solution to (1.2) with fo 6 L 1 ( M ) n 
there exists Cq> 0 depending on a, b, A,/o,£o and ||/o||i 1 (JVf) 


2 ), estimate 
such that 


(2.3) holds. Moreover, 


(3.1) 


sup| f{ft)\ < Cq 

t> 0 


Proof. We first apply Cauchy-Schwartz’s inequality to find 

(3.2) s '“i / s (X, / ) 1/2 - (X.= / ) 1/2 - 

for any / e P(R 2 ) n L 1 (w 2 ). Now, for f t a solution to ( |1.2[ ), we have 

iL f ‘ M = 

= } v? i l ~Ax-B e {/f t )r)ft- 
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Using the definition of A and B e , and then (3.21, we get 
d 


dt 


J ,Mdxdv < - J~ ( - 1 + ax 2 - bxv + v 2 (v - X)(v - 1 ) - ev 2 + xv + Iov^ft + £ ^/( ft ) 2 

f (v 4 + x 2 )f t + e f v 2 f t 

Jr 2 Jr 2 

f ft M dxdv, 

J R 2 


< A' , - K -2 

< A', - K 2 


where Ki and K 2 are generic constans depending only on a, 6, A, Jo and £q. Using Gronwall’s 
lemma we get (|2.3|) for some Co > 0. Finally, coming back to (|3.2|), we get 


□ 


\a?(ft )\ 2 < is, < 2 \\,\\l H m) ^ 2 max(C 0 , ||/ o || ii ( M )), 

which is nothing but S 

Lemma 3.2. For any ^ eR fixed, there exist some constants K\,K 2 > 0 depending on a , b , A, Jo, ^, k 
and £q such that 

(3-3) f Qeic/] f • sign(/) m < Jfi|/|i,i(R 2 ) - K 2 \\f\\ L i(m), V/eL^m). 

Proof. Since ^ e M. is now fixed, for simplicity of notation, we drop the dependence on this 
parameter. Using Kato’s inequality, sign (f)d 2 v f < d 2 v \f\, we have 


f Qe f ■ sign(/) 

Jr 2 


m < 


—k J ]/] (Ax + B e v - (1 + nv 2 ))m, 


thus 


f R2 Qe f ■ sign(/) m < - p(x , v) \f\ 5 


where p(x,v) is a polynomial on x and v with leading term v 4 + x 2 . Inequality (3.3) follows 
directly. □ 


Corollary 3.3. Estimate (2.4) holds. 


Proof. For , solution to (1.2), inequality (3.1) tells us that \^f(ft)\ < Cq. Moreover, since the 
mass is unitary for almost any t > 0, it holds 

d 


dt 




m < I\\ - K 2 


La 


t m , 


where Ki and K 2 depend only on a,b,X,I,£o and C 0 . Finally, integrating this last inequality, we 
get 

II ft Hi 1 (m) ^ max(Ci, |/o | z,i ( m )), Vt>0, 

for some positive constant C\ depending only on the parameters of the system, £q and C, □ 


Now we analyse the jU(m) and H 2 {m) norms of the solutions to (1.2), in particular, we prove a 
priori bounds (2.51 and (|2.6|). Since the equation is hypodissipative, we used the ideas of “twisted 


spaces” and the Nash-Villani’s technique (see e.g. m\) to control the L 2 (R 2 ) contributions in 
function of the L 1 (M 2 ) norm. 

Lemma 3.4. For 0 < k± < n 2 , let us consider two exponential weight functions m± and m 2 as 
defined in (£3. For any ^ € K. fixed, there exist Ki,K 2 > 0 and S e (0,1) constants such that 

(3.4) (Q e [^]/,/>^<^ 1 |/lli 2 ( R 2 ) -^ 2 ||/|^, V/ e U 1 (m 2 ), 

where (v)'H 1 the scalar product related to the Hilbert norm 

WfS : = \\f\\lHrn 2 ) + hS\dxf\\l2 {m2) + s 4/ ^d x f,d v f) L 2 (mi) + s\\d v f\\ 2 LHm2) . 
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Remark 3.5. It is worth emphasising that for S e (0,1) the norm H l is equivalent to the usual 
norm ofH 1 (m 2 ). Indeed, the choice of the exponents allows us to write 

e5/3 c 

c <5|l/llffi(m 2 ) ^ II/II L 2 (m 2 ) + 7j~) ll<7r/|lL 2 (m 2 ) + ^ ll^v/lL 2 (m 2 ) - ll/llw 1 ) 

for some c$ > 0. 

Proof. The proof is presented as follows: the first three steps deal with inequalities in L 2 for / 
and its derivatives, while the last one combines these inequalities to control the P 1 norm. Some 
long and tedious calculations are only outlined for the sake of clarity. In the following we denote 
by jfeo.jfei and k 2 some unspecified constants and drop the dependance on J?. 

Step 1. L 2 (m 2 ) norm. We start by noticing that 

(dvvfJ)L*(m 2 ) = - [ (dyffml + k 2 [ (l+ 2n 2 v 2 )f 2 m 2 2 . 

(d x (Af),f) L 2 (m2) = ^ f \d x A-A^h.y 2 m l = I f [a-2K 2 x(ax-bv)]f 2 ml, 

A JR* 777/2 ^ 

and similarly 

(d v (B e f),f) L 2 (m2) = \ f [3v 2 -2(1 + X)v +X-£-2k 2 v B e \,f 2 ml- 
Therefore, we get 

(3-5) ( Qef , /)i=(m 2 ) = ~ f R2 P( x ’ v )f 2 m 2 - ll^/li 2 (m2) , 

where p(x, u) is a polynomial in a; and v with leading term t; 4 + x 2 . In particular, there exist some 
positive constants k\ and k 2 such that 

(3-6) {Qef 1 /} L 2 (m 2 ) ^ hi ||/Hl 2 ( R 2 ) - k 2 \\f\\ i 2 ( M i/ 2 m2 ) - \d v f\\ LHm2 y 

Step 2. x-derivative bound. We have 

(d x (d x (Af)),d x f ) L , (ma) = i f hd x A-A^](d x f) 2 m 2 2 

v 7 2 JR 2 L 777-2 J 

= i r [3a-2^(ax-6u)](4/) 2 m^ 

and 

(d x (d v (B e f)),d x f) L 2 (ma) = /^[dvB&f + 9 x B e aj + B s d 2 xv f]d x fm 2 . 

Since = 1, and observing that 

f^[d v B e d x f + B e d 2 xv f]d x fm 2 = i £ 2 [d v B e - B e ^^]{d x f) 2 m 2 , 

we get 

{d x (d v (B s f)),d x f) L 2 (to2) < r |4/||^/|TO| + ' f \d v B e - B e d h rr ^](d x f) 2 ml. 

777-2 

Using that 

( d x d 2 v f,d x f) L2(m2) =- f R2 \d 2 xv f\ 2 m 2 2 + X - f R2 (d x f) 2 d 2 v m 2 2 . 

we finally obtain 

(3.7) {dx(Qef)id X f)L 2 (m 2 ) ^ fcl II^/IIl 2 ( R 2 ) - k 2 \\9 x f\\ L 2^ M i/2 m2 ^ 

-\\d X vf\\l 2 (m 2 ) + f R2 \dxf\\d v f\m 2 2 . 



10 


S. MISCHLER, C. QUININAO, J. TOUBOUL 


A similar calculation leads to 

(3-8) (( d v Q e f),d v f) L 2 (m2 ) < /ciR/lli 2 ( H 2 ) -k 2 \\d v f\\ 2 L2(M1/2m2) 

-\\dy V f\\l^m 2 ) + b f R2 \dxf\\d v f\m 2 2 

+ 2 ^ 2(1 +A) l~ v f 2 m 2 . 


Step 3. Cross product bound. The contribution of the cross product term is a little bit more 
delicate. We decompose it into five quantities and we study them separately: 


(0 x Qefi f)L 2 (m\) {^vQefi ^x/Jl 2 (mi) 

= [ 2 [( d x f){dl vv f ) + (d v f)(d%vvf)\™i 

JR Z 

+ f [d x Ad v f + d v Ad x f + Ad 2 x f] (d x f) ml 

Jr 2 l j 

+ f [d 2 vv Bf + 2 d v Bd v f + Bd 2 vv f]{d x f) m\ 

JR Z 

+ f R2 [2 d x Ad x f + Adl x f](d v f) ml 
+ [ [d v Bd x f + d x Bd v f + Bdl v f]{d v f)m\ =: £7". 

JR 2 i=l 

We start by handling the first term on the right hand side. Using integration by parts adequately, 
we get 

Tr = [ (d x f)(d v f)d 2 vv ml- 2 f(d 2 xv f)(d 2 J)m 2 . 

Jr 2 Jr 2 

Similarly, for the contributions involving A , we have 

T 2 = l f \d v A-A^](d x f) 2 m 2 + a f (d x f)(d v f)m 2 , 

2 Jr 2 l m\ -* Jr 2 


t 4 = 


f \d x A-A^](d x f)(d v f)m 2 + \[ d v [Am 2 ](d x f) 2 . 
Jr 2 l mf J 2 Jr 2 


Adding these last two expressions, it only remains 

, d x ml 


[ d v A(d x f) 2 m\+ f ha- A^ {d x f)(d v f)m\ < -b \\d x f\\ 2 L2(mi) + k 0 f \d x f\\d v f\M m\, 

JR? JR? TTl-\ JR 2 


for some constant ko > 0. 

For the contributions related to B E , involved in 73 and 7s, we have 

7-3 = - f2 Kl x(3v-l-X)f 2 ml + 2 f d v B £ (d x f)(d v f)m 2 + f B e {d 2 vv f){d x f)m\, 

and 

r 5 = d v B e (, d x f){d v f ) ml+ l - [< %B e - B E d ^](d v f) 2 ml, 

Finally, for the last contribution in 75, we have 


f B E (d 2 J)(d x f)m 2 < k 0 f (d 2 vv f){d x f)M 3 ' 2 ml 

JR 2 7 r 2 
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getting that there exists ko > 0 such that 

(3.9) ( dxQef , d v f) L 2 (mi) + (d v Q e f, d x f) L 2 (mi) 

<fc 0 f \d x f\ \d v f\Mm\ + ko f \d 2 xv f\ \d 2 vv f\ m 2 

-b Il^/I 2 i2(mi) + fco |0„ 2 „/l |3 X /| M 3 ' 2 m 2 
+fco f \d v f\ 2 M 2 m\ + k 0 f fMml 


Step 4- Conclusion. To get (3.4), we just put together (3.6), (3.7), (3.8) and (3.9) and we use 
Young’s inequality several times. Indeed, the scalar product (•,•)«! applied to any / e i7 1 (m 2 ) 
writes 


{Qef: f)'H x = (Qe/,/)L 2 (m 2 ) 

+8 3,2 {d x Q e f , d x f) L 2 (m2) + <5 ( d v Q e f,d v f) L 2 (m2) 

J 4 / 3 6 4 / 3 

+ — (d x Qefid v f) L 2 (mi) + —(d v QJ,d x f) L 2 {mi) . 

To give an idea of the method, we only explain how to get rid of a few terms. For example, for 
the positive contribution of ( |3.7[ ), it holds 

<5 3/2 M^/ll! 2 (R 2 ) + <5 3 ' 2 f \d x f\\d v f\ml < 11^/11^2^2) + S 7/4 R/|| 2 2 (ro2) + |a„/||| a(m2) , 

and for 5 > 0 small enough these terms are annihilated by the quantities 


-\\dvf\\l^ m2) -S 3/2 k 2 \\d x f\\ 2 LHM1/2 


8 4 ' 3 b 


m 2 ) 


\dxf\ 


L 2 {m 1) 


present in the right hand side of (3.6), (3.7) and (3.9). 


In (3.8), the only delicate contribution is 

8 3 ' 3 b 


Sb 


f \d x f\\d v f\ 


m 2 2 < 


\d x f\ 


8 1 Cb 


L 2 (m 2 ) 


l^/|lz, 2 (m 2 ); 


but the right hand sides of (3.6) and (3.7) include 


ll < 5t'/lli 2 (m 2 ) b 3/2 k 2 \\d x f\\ 2 L2(Ml/2m2) , 


and once again for 8 > 0 small the sum is nonpositive. 

The positive part of (|3.9[) is controlled using that < k 2 . Indeed, in that situation 


5 A ' 3 k 0 f \d 2 vv f\\d x f\M 3 / 2 ml < 5 4/3 - 1/4 fco||Si/||i 2(m2) +<5 4/3+1/4 fc 0 |S ;c /|| 2 i2(m2) , 

replacing, if necessary, ko by a larger constant. If 8 > 0 is small we get rid of these terms thanks 
to the presence of 

-8 3/2 k 2 \\d x f\\ 2 L2(M i /2m2) - 8 ||^/||| 2 (m2) , 

in (3.7) and (3.8). 

All remaining positive contributions can be handled in the same fashion leading to the conclusion 
that one can find K \, K 2 > 0 such that 


(Q e / J /)« 1 <A' 1 ||/iii 2(M2) -^ 2 |/i^ 1 . 

□ 


Corollary 3.6. Estimate (|2.5|) holds. 
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Proof. Nash’s inequality in the 2-dimensional case reads: there exists a constant C > 0, such that 
for any / € L 1 (K 2 ) n H 1 ' 

(3.10) 


C S' 

ll/lli 2 (R 2 ) - C'II/IIl 1 (R 2 )II- D x,z;/||l 2 (R 2 ) ^ A77 ll/lliqR 2 ) + ~^\\-^x,vf\\L 2 (R 2 )- 


26' "^ K ‘> 2 

Coming back to the previous lemma, using the equivalence of the norms 'H 1 and H 1 (rri 2 ), together 


with the fact that a solution f t to (1.2) is a probability measure, we get that, 


= (Qe[Sf t ]ft,ft)w * kl-kilftlw* 

for some Aq, k 2 > 0 constants. Finally, integrating in time, we get 

l/t||«i <max(C 2 , WfoWn 1 ), 

for some C 2 > 0 depending only on the parameters of the system and the initial condition. □ 
Let us notice that we can go a little further in the analysis of the regularity of the solutions 


of (1.2). Actually, we can expect that the norm H^(rn) is also bounded. Indeed, there exists ko > 0 
such that 


(3.11) 


J vv nL 2 


{9y V Q S f,d^ v f) L 2( ma ) 

= - f \dl vv f\ml + \ J \d 2 vv f \ 2 d 2 v 

! f (d v A)(d 2 vv f)(d 2 xv f) ml + i f \d 2 vv f\ 2 [d x A - A^} i 
■ f ( a v VV B)f(dl v f) ml+ 3 f (dl v B)(d v f)(dl v f) m 2 2 


+2 


1 

H— 
2 


f \dlf\ 2 [ 


5 d v B e -B t 


d v ml - 

\m\ 

ml - 

ml + j 

f I/I 


k o[f \ d vvf\ 2 + J \ d lvf\ 2m l+ f \f\ 2m l + f \ d vf\ 2 ml]. 


We can therefore state that 


Corollary 3.7. Estimate (2.6) holds. 


Proof. The proof follows the same idea already introduced in the proof of Corollary 3.6 We 
consider the norm 

WfWui : = ll/ll W 2 lOli 2 (m2) , 


and notice that (|3.6|), (|3.7|, (|3.8|) together with (|3.1l| imply that 

ftlm 


— II f 1,2 


<- ^ii/.ib 


+2« 2 t„[ f + f + f |/,| 2 ml+ f |a„/J 2 ra 2 ] 

< fcl - fc 2 ||/t||^2, 


for some k \, fc 2 > 0 depending on some <5 > 0 small and the parameters of the system. Inequality (2.6) 
follows. □ 


3.2. Entropy estimates and uniqueness of the solution. Now we focus our attention on the 


problem of uniqueness of the solutions to (|1.2|). First, we prove that solutions remain in the space 

:ive fun 

■ dxdv, 


of functions with finite entropy. To that aim, for any positive function /, we define 

\d v f(x,v)\ 2 


IM) ■- L 


fi.X,v) 

which is understood as a partial Fisher information. When the previous quantity is not well defined 
we use the convention /„(/) = +°°. Notice that in any case /„(•) > 0. Equipped with this definition 
we can state: 
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Lemma 3.8. For any /o 6 L l (M) n L 1 log L 1 nP(t 2 ) we denote by f t the associated solution to 
the FhN statistical equation (1.2) with initial condition / 0 . It holds 

(3.12) sup f* I v (fs)ds<C(T), 

te[0,T] “'0 

where C(T) depend on /o and the coefficients of the problem. 

Proof. It is well known that for functions with finite moments, the entropy can be bounded from 
below. Indeed, since 

ri logri > -r 2 + ri logr 2 , Vrr > 0, r 2 > 0, 
taking r\ = f(x,v) and r 2 = e~ M , it holds 


0 > / log/ > - e- M -fM, 


Jt?(ft)>- f e M - f f t M>-2ne 1 - max(C 0 , ||/o||ii(M))- 
J h 2 Jr 2 v ' 


implying that 


On the other hand, for any solution of ( |1.2[ ) with initial datum /o there exists a positive constant 
C', depending on the parameters of the system, £q and Cq, such that 

= f (1 +log( f t )) Q E [M ft 
= ~Wt)+ J (d x A + d v B e ( f fl )) f) 

^ -Iv(ft) + C\\ft\\ L i(M)- 

Let us fix T > 0 and take any t <T, thanks to estimate (2.3), we get that 

J4?(ft) < ~ I v (fs) ds + J'f’(fo) + CT max(C 0 , ||/o||li(m))- 

Since Jf? is bounded by below, we get that I v (ft) e L 1 ([0,T]). Moreover, taking the supremum 
on the last relationship, we get 


sup Jif(ft) < Jtf’(fo) + CT max(C Q , ||/o||li(m))- 

te[0,T] 


□ 


Corollary 3.9. For any two initial data /o, 5o € L 1 (M 2 )n L 1 logl/n P(R 2 ) the associated solutions 
f and g to the FhN statistical equation (0, satisfy 

SU P II ft ~ fl l t|li 1 (M) ^ C(T) |/o - go Hi 1 (M) i 

[0,T] 


for some positive C(T). In particular, equation (1.2) with initial datum in L l (M 2 ) n ZdlogL 1 n 
1)2 ) has, at most, one solution. 


Proof of Corollary \3.9\ We write 

dt(ft-gt) = Qe[/{ft)]{ft~gt)+e/{ft-g^dygt 

from which we deduce 

i [\ft-9t\M < K 1 [ \f t -g t \M + e\ c /(f t -gt)\[ \d v gt\M 

dt Jr 2 Jr 2 Jr 2 

< K x f R2 1 ft - g t \M + e\I(g t n t C {M2) \ft - g t |M, 

where A/ is the constant introduced in the proof of Lemma [3.1| Also, it is not hard to see that 
SU P IgtllzPCM 2 ) ^ IgolU^M 2 ) + 2(Ai + 1)Tmax(C 0 , IgolU^M))• 

te[0,T] 

The rest of the proof is a direct application of the time integrability of I v (gt) and Gronwall’s 
lemma. □ 
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Let us finish this section by giving some insights of the proofs of the existence of solutions and 
stationary solutions to equation ( 1 . 2 ) which are, however, classical. 

Proof of Theorem 2.2 Let us consider an exponential weight m and J! e L°°(R+) such that 


sup |^| < Co, 

t> o 


where C' 0 is given by (3.1). First, to avoid the non boundedness of the coefficients of the equation, 
let us fix R > 0 , and define a regular truncation function 


(3.13) 


Xr(x,v) =x(x/R,v/R), x 1 


')> l.B(0,l) ^ X ^ 1b(0,2)- 


Secondly, to avoid the intrinsic degenerate character of ( |1.2[ ), we fix some 1 > er > 0, and define the 
bilinear form 


^cr * {dvf-)dv9)L 2 (m)+{dvfi9XRTTl d v 7Tl 

{&xfi ^x9)L 2 (m) ® {&xf: 9 X.R ^x^R )L 2 (m) 


-^(f,9XR[9 x A- Am 2 d x m 2 ]) L 2 (m) 

-^{f,gXR [dyBe(J r t) - Be(ft) m~ 2 d v m 2 ]) L 2 (m) . 

This bilinear form is obviously well defined, a.e. t > 0, for any /,j t H 1 (m). Moreover, a a is 
continuous, 


| f: (?)| — Cr ||/|| H 1 (m) ll^/ll H 1 (m) ) 


for some positive constant Cr, and coercive. Indeed, we have from (3.6), that 


1 " ~ -112 , 17 11 n J-||2 


a<j(t;f,f ) > - \\dvf\\ L 2 (m) + - \\dxf\\ L 2(m) - M/llz, 2 (m)> 

for some k\ > 0 not depending on t, nor on R and nor on er. The J. L. Lions theorem [ 8 ] Theorem 
X.9] implies that for any /o e L 2 (m) there exists a unique 


/ 6 L 2 ((0, 00 );IL 1 (to)) n C([0, °°);L 2 (to)); e L 2 ((0, oo)-H\m)') 

at 

such that /( 0 ) = /o and 

(^ t f,g)L 2 ( m ) + a <7 (f(t),g) = 0, a.e. t > 0, V <7 e U 1 (to). 

We recall that /_ := min(/, 0) belongs to H l {m ), therefore we can use it as a test function to 
find that 

/o > 0 =>• f(t) > 0 , a.e. t > 0 . 

Let us now fix some T > 0. Using / itself as a test function, we get easily that 

WftWlRm) + f Q \\9vfs\\L 2 (rn)ds ^ efclT ||/o||| 2 ( m )> 

therefore, one can take the limits a -*■ 0 and R -+ 00 , to find that for any ip e C 1 ([0,T]; C 2 (K 2 )) 


f Ttft = f Tofo + f f [d t tp s + d 2 v tp s - Ad x tp s - B e (J r s )d v ip s ]f s ds, 0 < t < T, 

Jr 2 Jr 2 Jo Jr 2 l 

holds. Taking a well chosen sequence ip n -* Af 2 , we deduce that 

sup ||/t |z,i(M 2 ) < nrax(C", ||/o||li(m 2 )), 
tqo.r] 

for some positive constant C' that depends only on the parameters of the system. We also notice 
that, thanks to renormalisation concepts, we recover the inequality 

sup Jf’(/ t )+ f I v (fs)ds < J^(fo) + iLoTmax(C' 0 ,|/o| L i(M))- 
tqo.r] do 
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Let us take now f 0 e L 1 (M 2 ) n L l \ogL 1 nP(E 2 ), and a sequence {f n ,o} c L 2 (m) such that 
fn,o -*■ /o in L 1 (M). Moreover, let us assume that there is a positive constant C > 0 such that 
J^(f n ,o) ^ C, for any n € N. From the previous analysis we get a family {/„} e C((0,T); L 1 (M)) 
of functions related to the initial conditions {/ n ,o}- Using the Dunford-Pettis criterium we can 
pass to the limit in L 1 (M) finding a solution to the linear problem 


(3.14) 


dtf = d x {Af) + d v {B e {j? t )f) + d 2 vv f. 


that depends continuously to the initial datum (in the sense defined in Theorem 2.2). Moreover, 
from Corollary |3.9| we get that this solution is necessarily unique. 


Finally, we use again the ideas of Corollary 3.9 to find a solution to the NL equation (1.2). 
Indeed, it suffices to notice that the mapping 


I L°°([0,r]) — C([0,T];L 1 (M 2 )) 

1 S ~ /> 


with / solution of (3.14) for ^ given, is Lipschitz and contracting when T > 0 is small enough. □ 


Existence of stationary solutions will be shown as a result of an abstract version of the Brouwer 
fixed point theorem (a variant of US Theorem 1.2] and (TS]i: 

Theorem 3.10. Consider Z a convex and compact subset of a Banach space X and S(t ) a 
continuous semigroup on Z. Let us assume that Z is invariant under the action of S(t) (that is 
S(t)z € Z for any z 6 Z and t> 0). Then, there exists zq e Z which is stationary under the action 
of S(t), i.e, S(t)zo = Zq for any t> 0. 


We present the argument briefly in this section. Our aim is to find a fixed point for the nonlinear 
semigroup SQ e (t) related to equation (1.2). At this point we do not have any hint on the number 
of functions solving 

Qs[Jf]F = 0, 


and the nonlinearity could lead to the presence of more than one. However, in the disconnected 
regime e = 0 the nonlinearity disappears, and the multiplicity problem is no longer present. 


Proof of existence of stationary solutions to (2.7). Let us fix m an exponential weight and define 
for any t > 0 

S(t):X-+X with A = 7L 2 (m)nL 1 logL 1 nP(IR 2 ), 


such that S(t)fo is the solution to (1.2) given by Theorem 2.2 associated to the initial condition 
fo. Estimates (2.6) and (3.12) imply that S(t) is well defined. Moreover, the continuity of S in 
the Banach space L 1 (K 2 ) is direct from the definition of weak solutions, in particular, 

S(t)fo e C([0, oo); L 1 (M 2 )), 

with the topology of compact subsets in time. 

Finally, defining 


Z := Z(e) = {f e X such that (2.3) and (2.6) hold) c L i (K 2 ), 

which is invariant under S t for any t > 0 and convex. Moreover, the compactness of the inclusion 
Z c H 1 (m) L 1 (K 2 ) allows us to apply Theorem 3.10 and find the existence of a fixed point for 

S(t) and by consequence a stationary solution to (1.2). 

It is worth emphasising that the above proof show yet that the map e i 
in [0,oo), i.e., if £o > 0 is fixed, then 

G e e Z(e 0 ) for any e 6 (0,e o )- 


G e is locally bounded 


□ 
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4. The linearized equation 

The aim of the present section is to undercover the properties of the linearized operator associ¬ 
ated to Q e in the small connectivity case using what we call a splitting method. To illustrate the 
ideas we use, let us assume that an operator A on a Banach space X can be written as 

A = A + B, 

where A is much more regular than B, and B has some dissipative property. If B has a good 
localisation of its spectrum, under some reasonable hypotheses on A, we expect £(A) to be close 
to £(£). 

This is nothing but the Weyl’s abstract theorem (and/or the generalisation of the Krein-Rutman 
theorem) from Mischler and Scher [2Bj, that we recall here: 

Theorem 4.1. We consider a semigroup generator A. on a “Banach lattice of functions” X, and 
we assume that 

(1) there exists some a* € R and two operators A,B € 'if(X), such that A-A + B and 

(a) for any a > a* ,£ > 0 , there exists a constant C a / > 0 such that 

Vt > 0, | S B * CASB) ( *°(i)ll»( a) < C a/ e at . 

(b) A is bounded, and there exists an integer n> 1 such that for any a> a*, there exists 
a constant C a ^ n > 0 such that 

V i > 0, \\(AS B )(* n \t)U (x , Y) < C a , n e at , 

with Y c D( A) and Y c X with compact embedding; 

(2) for A* the dual operator of A defined in X', there exists /3 > a* and if e D(A*) n X' + \ {0} 
such that 

A*if > 13if ; 

(3) S\ satisfies Kato’s inequalities, i.e, 

V/€D( A), A6{f)>e'U)Af, 

holds for 9(s ) = |s| or 9(s ) = s+. 

(4) -A satisfies a strong maximum principle: for any given f and 7 e R, there holds, 

I/I € D( A) \ {0} and (-A + 7 )|/| > 0 imply f > 0 or f < 0. 

Defining 

A := s(A) = sup{i?e(/) : £ e £(A)}, 
there exists 0 </„£ D(A) and 0 < (f e D(A*) such that 

A / 00 = A/ 00 , A*(f = \cf. 

Moreover, there is some a € (a*, A) and C > 0 such that for any /o e X 

||S A (t)/o - eM (fo,f>)U\\x < Ce m l/o - (/o,0)/oo|Ly. 


From Theorem 2.3 we know that for any value of e there exists at least one G e non zero stationary 
solution of the FhN kinetic equation (1.2). The linearized equation, on the variation h := / - G e , 
induces the linearized operator 


&eh = Qe{jf{G e ))h + s/{h)d v G e . 

Moreover, let us recall that in Section [3] we proved that 

(0e[/(G E )]/,/)L2(m) ^ ^l|/IU 2 (K 2 ) “ -^21| /|U 2 (m) ) 

if we could make K± = 0, then the operator Q e together with Jz? £ would be dissipative. Since it is 
not the case, let us fix a constant N > 0 and define 


(4.1) 


B e A£ e - A, where A = N xr( x i v )\ 
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with xr given by (3.13). We remark that A e 3§(H 2 (m)), and that Af vanishes outside a ball of 
radius 2 R for any / e H 2 {m). 


4.1. Properties of A and B e . We now precise the dissipative properties of «5? £ . In particular, 
we present two lemmas dealing with the hypodissipativity and regularisation properties of the 
sppliting A and B e . We use some ideas developed in [27i ns] and [26] . 

Lemma 4.2. For any exponential weight to, there exist some constants TV, R > 0 such that ( B e + 1) 
is hypodissipative in H 2 (m). 

Proof. From the characterisation of hypodissipativity given in Section [2j it suffices to show that 
there exists a constant C > 0 such that 

l%(OIU(flJ(m)) < Ce *, t > 0, 

or simply, to show that for any h 6 H 2 (m), it holds 

(4-2) (S E /l, <— | ft || H^(m) > 

for some norm |) • equivalent to the usual norm |[ • ||# 2 ( m ). 

Let us recall that the operator B e writes 

B e = JS? E - A = {QA/G e }-N X R)h + e/{h)d v G e , 

and since e K. is a real constant, we can use all a priori estimates on Q e directly. As usual, 
when no confusion is possible, we drop the dependence on a f e . Three steps complete the proof: 


Step 1. Dissipativity in L 2 (m). Let us notice that for any h € L 2 (m) we have 

\f{h)\<C\\h\\ mm) , 
for some constant C > 0. It follows that 

f{h) [ (d v G e )hm 2 < \/{h)\\\d v G s \\ L2(m) \\h\\ LHm) < C \\d v G s \\ LHm) f h 


2 2 
m . 


Thus, coming back to (3.5), we find that for TV and R large enough one can assume k\ - 
getting 

(4-3) (B e h,h) R 2(m) ^ ~ll^llz, 2 (m) ~ ^2 ||h|| i2 ( M i/2 m ) - ll^u^ll L 2 (m) > 

as a consequence, {B e + 1) is dissipative in L 2 {m). 


Step 2. Bounds on the derivatives of B e . For the ai-derivative we see that there exists some constant 
C' depending on X R and its derivatives, such that 

-N(d x ( XR h),d x h) L 2 {m) < C'\\h\\l 2(m) - N\\(d x h)^xE\\ 2 L2 {m) . 

On the other hand, thanks to Young’s inequality, we get 

J?{h) f {dl v G e ){d x h)m 2 = -f(ii) f d x G £ [d 2 x h + 2 kv d x h] rn 2 

< /{hf ||0„G e ||£ a(m) + \\\d 2 xv h\\ 2 mm) + \\V2Kvd x h\\l H my 

These two inequalities, together with ( |37t , imply that for TV and R large enough 
(d x (B £ h),d x h) LHm) < -\\d x h\\ 2 LHm) - \\\dl v h\l^m) + + L 


Proceeding similarly with the u-derivative we get 

/(h) f (d 2 vv G E ){d v h)m 2 = \^(h)\\\d 2 vv G4 L 2 (m) \\d v h\\L^ m) 

^ 2 II^TOGe||i 2 (m)(C 2 ||^ll 2 (m) + ll^^lli 2 (m))> 
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then, coming back to ( |3.8| ), we find N, R > 0 such that 

(d v (B s h),d v h) L 2 ( m) < -\\d v h\\l H m) - \\d 2 vM 2 L H m) + C' \\h\\ 2 L2{m) + jf a \d x h\ \d v h\ 


m 2 . 


Finally, for the second v- derivative we find C' such that 

-N(dl( X Rh),d 2 vv h) LHm) < -N [ XR (d 2 vv h) 2 m 2 + C' f (d v h) 2 m 2 + C'f h\d 2 v h\m 2 , 
and for any e > 0 

J R2 (9v VV G e )(d 2 v h) m 2 < ^— + e(\\$ v Ge\\i 2 (m)\\dvvvH 2 L 2 ( m ) + 

+ \\^vv^e\\L 2 (m) \\^ KV i^vv^ 1 ) lli 2 ( m ) ) • 

If e > 0 is small and N,R large enough, we obtain as an application of (3.11|, that there is a 
constant C' > 0 such that 

{dvv(B e h),d 2 v h)L2( m ) < -\\dv V h\\ 2 L2(m) + ^[ll^lli^m) + \\dvhf LHm) + \d 2 v h\ 2 L 2^ m ) + 


Step 3. Equivalent norm and conclusion. Let <5 > 0 and hi,h 2 e H 2 (m ), we can define the bilinear 
product 

{hi, h 2 )H%(m) := ( ^1, h 2 ) L 2 {m) + ^{d x hi, d x h 2 ) L 2 (m) + S{d v hi,d v h 2 ) L 2 ( m } + ^ {9 xv hi, d 2 v h 2 ) L 2 (m) ■ 

and the relative norm 

Ill'll ffj(m) := II ^ II L 2 (m) + 8 \Dx,vh\ L 2( m } + S \\d vv h\\ L 2( m y 

Choosing <5 > 0 small enough we conclude that for any a e (0,1] one find S a such that 

{B e h, h)H 2 (m) - _Q! Ill'll ~H 2 {m)' 

Since the norm related to H 2 (m) is equivalent to the usual norm in H 2 {m), we can conclude that 
(. B e + 1) is hypodissipative in H 2 (m). □ 


Lemma 4.3. There are positive constants N,R large enough and some Cjs e > 0, such that the 
semigroup Ss e satisfies 

\SB e {t)h\ H 2( mi } <CB e t 9 / 2 |/l| L 2( m2 ), Vie (0,1]. 

As a consequence, for any a > -1, and any exponential weight m, there exists n > 1 and (7 n>e such 
that of any t > 0 it holds 

(4-4) \\(AS Be )(* n \t)h\\ H 2 (m) < C n , e e at \\h\\ L 2 {m) . 

Proof. We split the proof in three steps, in the first one we refine the previous estimates on the 
norm of the semigroup associated to the operator B e , in the second one we use Hormander-Herau 
technique (see e.g. [50]) to get the first inequality, and finally we prove (4.4). 


Step 1. Sharper estimates on B e . We denote for K > 0 a generic constant. From the proof of the 
previous Lemma, we know that there are N, R large enough such that for any h € D(B e ), it holds 

{B e h,h) L 2( m2 ) < -K\\h\\ 2 L 2( m2 ) - ||^ft.|| 2 2( rn2 ) 


1 ’ II 2 1 II »2 > || 2 , jy || i || 2 [ 


{d x B e h,d x h) L 2( mi ) ^ 9 ll^x^|li, 2 (mi) 9 ll^xi;^lli 2 (m 1 ) + -^I^Ib 2 (toi) + 2 ^JI^^Hi 2 (TOi) 


(d v B e h,d v h) L 2 ( mi ) < ~\dt v h\ 2 L 2( mi ) + K\\h\\ 2 L 2^ mi ) + — \\d v h\\ 2 L2 ^ mi ^ + St\d x h\ 2 L2 ^ mi ) 


St 

1 


(d vv Beh,d vv h) L 2 ( mi ) < K\h\ L 2( mi -) + K\d v h\ L 2^ mi ) + ^\\9 xv h\\ L 2^ mi y 
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We also notice for any 5,t e (0,1) it holds 


((Qe + {^v(Qe ^)^x^)L 2 (mi) — ^ ll^c^'iz/ 2 (mi) 

+ K ll^lll,2( m2 ) + — ||3„/l||i2( m2 ) + JpjTo II^TO^Ili 2 (m 1 ) + ^^ 1/10 ||^/i||^2( mi ), 

(f{h)d 2 xv G ei d v h) LHmi) + {j?(h)d 2 vv G s ,d x h) LHmi) < l9xGe l LHmi) [2S(h ) 2 


I II/} 2 ftll 2 I II/} /j || 2 1 | i 5 ™ Ge i i2 ( m i) rc/C^) 2 ,r||a u II2 ] 


yielding to 


fo A' 

(d x B e h,d v h) l 2( m ) + (d v B e h,d x h) l 2( m ) < --||9 x /i|| i 2( mi ) + ^ ||/i|lx,2( m2 ) + '^ll l %^llL 2 (m2) 

A 


K 

tS 


^1710 R 2 „fc|li2 (mi) + A^ 1/10 ||/} 2 >|| 2 2 (mi) . 


Step 2. Hormander-Herau technique. For a given h e H 2 (mi) n A 2 (m 2 ) we denote /i* := Sg E (t)/i, 
and define T by 

•TXM) := \\h\\l 2(m2) + Cl t 3 \\d x h\\l 2(mi) + c 2 t\\d v hf L 2 (mi) + c 3 t 2 (d x h,d v h) L 2 (mi) + c 4 t 4 \\d 2 v h\\ 2 L 2 (mi) , 

which, for well chosen parameters, is decreasing. Indeed, thanks to the inequalities found in the 
first step, we have 

j t HtM < 177, 


with 


Tl 

t 2 

r 3 


= A' 


f \ - 2m\ + 2(cif 3 + c 2 t + c 4 f 4 )m? + ]/i?, 

Jr 2 l o 

f [(3ci + 2c 2 i5 - ^c 3 + 2c 3 5)t 2 - cit 3 ](a a ,/i t ) 2 m 2 , 
w/i? 2 4 

[ [ - 2m 2 + c 2 m 2 + + ^m 2 + 2c^Km\ + C -^ m 2 2 ](d v h t ) 2 , 

Jr 2 l odd d 

T 4 = f R J 3 [- c i + J +c 3 K6 1/w ](dl v h t ) 2 ml, 

Choosing 

c 4 = (5 2 , c 2 = 6 4/3 c 3 = S 312 and c 4 = <5 4 , 
we get that for Se (0,1] small enough, it holds 

j/(t A) < o. 

for any t e (0,1]. Since 0 < c 4 < c\ < c 3 < c 2 and c 4 c 2 > c§, we finally get that 

c±t 9 l 2 [\d XtV ht\ 2 L 2( mi ) + \\d% v ht\\ 2 L 2( mi ' ) ') < F(t,ht) < F(0,h o ) = ||/io||z,2( m2 ). 

Step 3. Proof of inequality ( |4.4[ ). From the definition of A we notice that 

\\ASB e (t)h\\ H 2 (m ) < C"t“ 9/2 e _t \\h\\ L 2 (m) , Vie (0,1], 

for some constant C'. It is important to remark that since A lies in a com pact , we do not need 
anymore two different weights m 4 and ?n 2 . Therefore, we apply Proposition 2.5 with X = L 2 (m), 
Y = H 2 (m), 0 = 9/2 and a* = -1 to get (4.41. □ 
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4.2. Spectral analysis on the linear operator in the disconnected case. We consider in 
this section the disconnected case e - 0. The corresponding FhN kinetic equation is linear and 
writes 

d t g = d x (Ag) + d v (B 0 g ) + d 2 lv g 
Bq = v (v - A) (v - 1) + x, 


Theorem 2.3 states that there exists at least one function Go e Pn Hf(m) which is a solution to 
the associated (linear) stationary problem 

J%G 0 = d x (AGo) + d v (B 0 G 0 ) + d 2 vv G 0 = 0. 

Since the operator now enjoys a positive structure (it generates a positive semigroup S^ 0 ), we 
can perform a more accurate analysis. Indeed, we can apply the the abstract Krein-Rutman 
theorem |4.1| previously stated. 


Proof of the stability around e = 0 in Theorem\2.3\ Let us assume for a first moment that hypothe¬ 


ses of the abstract Theorem 4.1 hold for Jz?o with a* - -1. We easily remark that 

A = 0, /oo = Gq 0=1, 

therefore, there exists a € (-1,0) such that 

£(.5?o)nA a = {0}, 

and 

^ /o € L 2 {m) , Vi > 0 (j Sse 0 {t)fo ~ (/o)Go||L 2 (m) ^ G e0t \\fo ~ (/o}Go||L 2 ( m )- 

Now, for e > 0, we consider G e such that 

Qe[J r G s ]Ge = 0, 

then, it holds 

j f {G E - Go) + J? 0 (G e -G 0 ) = h, h = sd v ((v - f(G e ))G E ), 

and, thanks to Duhamel’s formula, we get that 

II G e - Go\L 2 {m) ^ I Sp? 0 {t){G e - Go)||l 2 (to) + \\Spe 0 {t - s)/i||L 2 (m) ds. 

But G e - Gq and h have zero mean, then 

II Ge - Go||i 2 ( m ) < GflGe - Go||L 2 (m) eQt + e rrj II G e ||jyi(jvf i/2 TO ) (1 - e Qt ). 

Letting t ->■ oo we conclude that there exists C & > 0 such that 

||G e - G 0 ||i2 (m) <eG a ||G e || 

Finally, thanks to Corollary |3.6[ we have 

0 = (Q s [f Ge \G E ,G e ) H , < Ki - K 2 \\GeW^i < K^csKiWGeWli^, 

for any exponential weight m 2 . If /-c 2 > k, we have then 

ll^lffRM^m) - C K , K2 \\G £ \\ 2 Hl(m2) < C^KjcgKi, 

and in the small connectivity regime £ e (0,£o), constants I\\ and K 2 do not depend on e. Defining 
g(e) = £G a G KiK2 /\i/ciA' 2 we get the stability part of Theorem 2.3 


It only remains to verify that the requirement of Theorem |4.1| are fulfilled for i n the Banach 

lattice X = L 2 (m). 
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(1) (a) the splitting (4.1) has the nice structure. Indeed, the Lemma 4.2 implies that Bq + 1 

is hypodissipative in L 2 (m), therefore 

\\SB 0 (t)\\ss(L 2 (m)) £ Ce~\ V t > 0, 
i.e., it suffices to take a* = -1. 

(b) if Y = H 2 (m) and X = L 2 {m), the desired inequality is consequence of Lemma 

(2) The requirement is obtained for /3 = 0 and i[> = 1. Indeed, in that case 


4.3 


^ U 


(3) A side consequence of (3.6) is the positivity of the semigroup: 

/o>0 => Sjf 0 fo(t) > 0, Vf>0. 


Moreover, using that L 2 {m) is also a Hilbert space, we deduce the Kato’s inequalities. 

(4) The strict positivity (or strong maximum principle) is a straightforward consequence of 
Theorem |B.l | in Appendix |Bj 

□ 


Let us finish this section by summarizing the properties of the spectrum of A? 0 in the Banach 
space L 2 {m) and by a useful result on the regularisation properties of 1Z^ 0 (z). 


Proposition 4.4. 

(i) There exists a < 0 such that the spectrum £(2z?o) o/2zf)b in L 2 {m) writes 

E(.i? 0 )nA a = {0}, 

and 0 is simple. 

(ii) For any a > a, there exists a constant C#i > 0 depending on (a - a), such that 


Pl% 0 (z)\ £?(l 2 ( m),Hl(m )) — c Hi (i + M' 1 ), 

Proof. It only remains to prove (ii). Let us consider z e A a 


V z e C \ {0}, Re{z) > a. 

\ {0}, and take f,g e L 2 (m) such that 


(-2o ~z)f = g. 

Thanks to Lemma |4.2| and the definition of A, we get 

(Re(z) - d)|/||| 2 ( m ) + \d v f\ 2 L 2 ^ m ) < |s , IU 2 (m) l|/lli 2 (m) + AT ||/||| 2 ( m ). 

Moreover, (i) tells us that 0 is an isolated simple eigenvalue for Afo in L 2 (m ), then TZ^ 0 (z) writes 
as the Laurent series (see for example [23J Section 3.5]) 

oo 

n^ 0 (z)= ^ z k C k , C k € £3{L 2 {m)), 

k =-1 

which on a small disc around 0 converges. Thus, there is some C° > 0 such that \1Zse 0 ( 2 ) ||^(L 2 (m)) ^ 
C° |z| _1 for any z 6 A a , z t 0. Finally, we notice that 

min(l,a-d)||/|| H i (m) < (1 + 7VC'°|z|“ 1 ) ||ff||z, 2 (m 2 ) , 

therefore, it suffices to take Ch 1 = 1 + /NC° min(l, a - a), with N large enough. □ 


5. Stability of the stationary solution in the small connectivity regime 

Now, we establish the exponential convergence of the nonlinear equation. To that aim, we first 
notice that, in the small connectivity regime, the linear operator T£ e inherits (in a sense that we 
precise later on) the stability properties of Jz?o- 
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5.1. Uniqueness of the stationary solution in the weak connectivity regime. As a first 
step in the proof of Theorem |2.4[ we need a uniqueness condition that, for instance, can be settled 
as a consequence of the following estimate: 


Lemma 5.1. There exists a constant C\> such that for any g e A 2 (to), (g) = 0 and for the solution 
f e Lr{m) to the linear equation f = g there holds 

(5-1) ||/||v := II / II L 2 (Mm) + |V„/|| i 2 ( M i/ 2 m ) < Cv |$ II L 2 (m) ■ 


Proof. We easily compute 

f (J^f )fMrn 2 = - f p(x,v)f 2 m 2 - f ( d v f) 2 Mm 2 , 

J r 2 Jr 2 Jr 2 

for some p(x,v) polynomial in x and v with leading term v e + x 4 . Therefore, there exists some 
constants K\ > 0 and 0 < K 2 < 1, such that 

[ (J? 0 f)fMm 2 < K\ f f 2 m 2 -K 2 f f 2 M 2 m 2 -K 2 f ( d v f) 2 Mm 2 . 

Jr 2 Jr 2 Jr 2 Jr 2 


The invertibility of Jz?o in L 2 (m) for zero mean functions, writes 

■^of = 9 £ L 2 (m), (g) = 0 => \f\ L 2 (m) ^C a \g\ L 2 (m) , 

with C & given in the proof of the stability part of Theorem |2.3[ As a consequence, for any / and 
g as in the statement of the lemma, we have 


f f 2 M 2 m 2 + f (4/)W<-} f gfMm 2 + I p- f f 

Jr 2 Jr 2 A 2 Jr 2 A 2 Jr 2 


2_2 

m 


< - 
“ 2 


- f f 
2 Jm? j 


2 M 2 m 2 + 


2 K* 


[ 

Jr 2 


2 2 

gm + 


Kn 


f 

Jr 2 


2 _2 
g m , 


from which (5.1) immediately follows. 


□ 


Corollary 5.2. There exists e± e (0,£o) such that in the small connectivity regime e e (0,£i) the 
stationary solution is unique. 


Proof. We write 


G e -F e 


(5.2) 


e^ 0 -'[O v ((v - J(F e ))F £ -(v- e /(G E ))G e ]\ 
e^ l [d v ({v - f{F e )){F e - G e ) + (J(F e ) - S{G e ))G e )\ 


As a consequence, using the invertibility property of for zero mean functions, and the uniform 
bound (2.5) on G e , F e , we get 


ll^-G £ ||v 


< e c a ||5„((u - /(F e ))(F e - G e ) + (j?(F e ) -/{G e ))G e )|| L2(m) 

< eC\\F s -G e \\ v , 


for some C depending on the parameters of the system and £q. The previous relationship implies, 
in particular, that ||F £ - G e ||y = 0 for e < e\ = 1/C. □ 


5.2. Study of the Spectrum and Semigroup for the Linear Problem. We now turn into a 
generalisation of Proposition |4.4| in the case e > 0 small. Since the positivity of the operator is lost, 
Krein-Rutman theory does not apply anymore, however we can prove the following result based 
on a perturbation argument 

Theorem 5.3. Let us fix a e (a,0). Then there exists e 2 e (0,£i) such that for any e e [0,e 2 ], 
there hold 

(i) The spectrum S(-Sf £ ) of J£ e in L 2 (m) writes 

E(if £ )uA Q = { Me }, 
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where p £ is a eigenvalue simple. Moreover, since J? e remains in divergence form, we still 
have 

Jzf £ * 1 = 0 

and then p £ = 0. 

(ii) The linear semigroup S^ e (t) associated to T£ e in L 2 (m ) writes 

s^(t) = ^ t n e + R e (t), 

where II e is the projection on the eigenspace associated to p e and where R £ (t) is a semigroup 
which satisfies 

\\Re(t)\\ag(L2(m)) ^ ^i 

for some positive constant Cse ei independent of e. 


To enlighten the key points of the proof we present it in three steps: accurate preliminaries, 
geometry of the spectrum of the linear operator in the small connectivity regime and sharp study 
of the spectrum close to 0: 

Step 1. Accurate preliminaries: Let us introduce the operator 

P e = J? £ - i ? 0 = -ed v ((v - /(G e ))-) + £/{■) d v G e . 

Our aim is to estimate the convergence to 0 of this operator in a suitable norm. We notice that, 
for two exponential weights mi,m 2 as in ( |2.1[ ) with < K 2 , it holds 

\\Peh\ 2 L , {mi) < Ce 2 f K 2 {h 2 + v 2 \d v h\ 2 )m 2 1 + Ce 2 J?(h ) 2 

* C£2 (\\h\\ 2 L Hmi ) + \\dvh\\ 2 L 2 (m2) ), 

where C depends only on the parameters of the system and, in the small connectivity regime, on 
£ 1 . Therefore, there exists Cp ei >0 such that 

\\ P eh\ |l 2 ( mi ) < C Pei £\\h\\ ■ 


Step 2. Geometry of the spectrum of T£ e . 

Lemma 5.4. For any z e A a , z + 0 let us define K e (z) by 

K £ (z) = - P e K& 0 (z) AK Be (z). 

Then, there exists 772 (e)-*■ 0, such that 

£-*■0 

Vz«!J s := A a \ B(0, 772 (e)), \\K e (z)\\gg { L 2 {m)) <T 7 2 (£)(1 + 772 (e)). 

Moreover, there exists £2 € (0,£i] such that for any e € [0, £ 2 ] we have 

(1) I + K e (z) is invertible for any z e fl £ 

(2) T£ e - z is also invertible for any z efl £ and 

Vzen £ , 7 Zj? e (z) = U £ (z){l + Ji e (z)) _1 

where 

u £ {z) = n Be (z) - Hj? 0 (z) An Be (z). 

We thus deduce that 

E(^)nA a cB(0, 772(e)). 
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Pro of. W e define mi and m 2 two exponential weights with m 1 = m. From Lemma |4.2[ Proposi¬ 


tion 


4.4 


and the Step 1 we get that for any z € any /i € A 2 (m) 

A' e (z)/i|| L 2( m) < eC Pei [|^ 0 (2;)^l^ Bc (z)/i]| ff i (m2) 

< eCp ei Cpi(l + \z\~ 1 )\\AU Be (z)h\\ L 2 ( rn2 > ) 

< eCp £i (7pi(l + M _1 ) C ei \\h\\ L 2(m), 


where C ei is an upper bound of \\ATZs e \\ss(L 2 (m),L 2 (m 2 )) and do not depend on e. Defining 

772 (e) := (e6V £] C Hl C ei ) [ l\ 

it holds 

\\Ke(z)\\^ L 2 {rn)) < ? 7 2 (e) 2 (l + 772 (e)" 1 ) = 772 (e)(1 + 772 (e)), Vz e 
therefore, fixing s 2 > 0 such that 


772 (e) < 1/2, Ve € (0,£ 2 ], 
we obtain the invertibility of I + AT £ (z). 

Finally, for any z e fl s : 

(Jz^ - z)U E (z) = I + AT e (z), 

then there exists a right inverse of 2z? e - z. The rest of the proof is similar to the proof of [.'161 
Lemma 2.16]. □ 


Step 3. Sharp study of spectrum close to 0. 

Let us fix r e (0,-a] and choose any e r e [ 0 , 62 ] such that ri 2 (£r) < r in such a way that 
£(Jz? e ) n A a c B(0,r) for any e 6 [0,e r ]. We may define the spectral projection operator 


IF :=• 


9 — f Kx> e (z')dz'. 

Z7TZ J\z'\=r 


We have then the 


Lemma 5.5. The operator Tl e is well defined and bounded in L 2 (m). Moreover, for any £ e [ 0,e r ], 
it holds 

IIn e - n 0 ||^ (L 2 (m)) < 773 (e), 

for some 773(e) -» 0. 

£—►0 

Proof. Let us notice that 

n ° = -' R '-Z‘oA'R. Bo ( z '))dz' = J^n^ 0 AlZ Bo {z)dz' 

and 

IF = -2^X,| =r (^B.(z')-^o A K Be (/))(/ + K s (z'))- 1 dz' 

= / *«.(*') Ke(z')(I + AF(F))- 1 dz' 

+ ^ [ Hjz 0 A-Rs e {z'){I + K e {z')y x dz'. 

ZlTZ J\z'\=r 

Then, we deduce that 

IF-IF = f ft Be (z')A7(z')(/ + A7(z , ))- 1 dz' 

Zltl J\z'\=r 

+ ^~ [ t Kjz 0 A{K B Sz')-n Ba {z'))dz' 

Zm J\z'\=r 

+ [ Hx 0 An Be {z'){I-{I + K e {z'))- x )dz', 

Zm J\z'\=r 




KINETIC FITZHUGH-NAGUMO SYSTEM 


25 


here, the first and third terms are going to 0 because of the upper bounds of K e (z). For the second 
term, it suffices to notice that 


n Bc {z')-n B o{z') = n Bo (z')(B £ -B 0 )n Bs (z'), 

and use that (B s - Bo) = P £ . 

To conclude the proof we recall the following lemma from [23] paragraph 1.4.6] 


□ 


Lemma 5.6. Let X be a Banach space and P , Q two projectors in such that \P-Q\gg(x) < 1. 

Then the ranges of P and Q are isomorphic. In particular, dim(R(P)) = dim(R(Q)). 

Provided with this lemma and fixing s' such that 773 (e') < 1, we get the 

Corollary 5.7. There exists s' > 0 such that for any £ e [0,e / ] there holds 

£(Jzf e ) n A q = {/ u £ } and the eigenspace associated to n E is 1-dimensional. 

5.3. Exponential stability of the NL equation. In the small connectivity regime e e (0,e'), 
let us consider the variation h := f e - G e , with f e the solution to (1.2) and G e the unique solution 
to (2.7) given by Theorem 2.3 By definition, h satisfies the evolution PDE: 

d t h = Jz? 0 h-£d v (vh) +£ c /(f E )d v f e -£ i /(G E )d v G E = Jz? e /i + £ c /(h)d v h, 

moreover, the nonlinear part is such that 

L 2 (m) — Ill'll L 2 (m) 11^^11 L 2 (m) 

for some positive constant C. 


Proof of Theorem \2-4\ Let us first notice that, thanks to inequality (2.5) and the definition of 
we have that 

||£ c /(/i)<9„/i|| L 2 (m) < C , ArLe|/i| L 2 (m) , V/i 0 e 

where 

Cnl = Q 1 max(C2, ||/io||ffi( m ))- 
On the other hand, Duhamel’s formula reads 

h = Ss? e (t)ho + J S& e (t ~ s)(e /(h)d v h) ds, 

then, we have that 

rt 

— II ^ IIL 2 (to) — ll^Sfe || L 2 (m) + ||*^^f e (^ — || L 2 (m) 

^ e at \\h 0 \\ L 2( m ) + Cj? ei C NL £ e“ (t “ s) ||/i| i 2( m) ds 

= e at u( 0 ) + C^ C Nh £ £ e a ^- s \i( s )ds. 

In particular, 

u(t) < G^ 1 u{Q)e {a+C ^ CNL£)t , 

Summarising, it suffices to define i]*(e) := C 2 /\/£ to get that for any / 0 such that 

l/o “ Ge||ffi(m) ^ V*( £ ), 

it holds 

\\fe(t)-G £ \\ L 2 {m) < C^Jf 0 -G £ \\ L ^ m) e aH , 

with 

a = a + Cjz ei Cg 1 C 2 '/eP < 0 , 

if £* is small enough. □ 
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6. Open problems beyond the weak coupling regime 

In the weak coupling regime, we have demonstrated that existence and uniqueness of solutions 
persist. In that regime, noise overcomes nonlinear effects and the system is mixing: one finds a 
unique distribution with an everywhere strictly positive density. As coupling increases, highly non¬ 
trivial phenomena may emerge as nonlinear effects of the McKean-Vlasov equation. For instance, 
it is likely that in another asymptotic regime in which coupling is non-trivial and noise goes to 
zero, Dirac-delta distributed solutions shall emerge (in which all neurons are synchronized and 
their voltage and adaptation variable are equal to one of the stable fixed point of the deterministic 
Fitzhugh-Nagumo ODE). 

Here, we numerically explore the dynamics of the Fitzhugh-Nagumo McKean-Vlasov equation 
using a Monte-Carlo algorithm. We observe that complex phenomena occur as the coupling is 
varied. That numerical evidence tends to show that several additional equilibria may emerge, the 
stability of stationary solutions may change as a function of connectivity levels, and attractive pe¬ 
riodic solution in time may emerge. These regimes are particularly interesting from the application 
viewpoint: indeed, among important collective effects in biology, from large networks often emerge 
bistable high-state of down-states (characterized by high or low firing rates), and even oscillations. 
These two phenomena are particularly important in developing and storing memories, and this 
occurs by slowly reinforcing connections \2‘2] . Interestingly, these two types of behaviors emerge 
naturally in the FhN McKean Vlasov equation beyond weak coupling. For instance, for fixed 
a = 0.5, we present the solutions of the particle system varying the connectivity weight beyond 
small values, both in the bistable case (in which the FhN model presents two stable attractors) and 
the excitable regime, the most relevant for biological applications, characterized by a single stable 
equilibrium and a manifold separating those trajectories doing large excursions (spikes) from those 
returning to the resting state directly. In both cases, we observe (i) that the unique stationary 
solution is not centered close from a fixed point of the dynamical system: neurons intermittently 
fire in an asynchronous manner for small coupling. As coupling increases, a periodic attractive so¬ 
lution emerges, before the appearance of distinct stationary solutions (two in the bistable case, one 
in the excitable case). These phenomena are depicted in Fig. [l] Proving, for larger coupling, the 
existence and stability of a periodic solution or distinct and multiple stationary solutions constitute 
exciting perspectives of this work. 

These phenomena are actually conjectured to be generic in coupled excitable systems subject 
to noise. 


Appendix A. Mean-Field limit for Fitzhugh-Nagumo neurons 


Let us start by a well known result with is a simple application of global existence and path wise 
uniqueness for system of SDE, see ni Chapter 5, Theorems 3.7 and 3.11] for example. Consider 
the particle system for 1 < i < N: 


(A.l) 


T JV 

dv\ = {v\ ( v\ - A) (1 - v\) - x\ + J 0 ) dt + — £ (vl - v° t ) dt + dW t l 

3 = 1 

dx\ = (■ -ax\ + bv\)dt, 


with initial data (Vq, Vf) for 1 < i < N distributed according to /o e P 2 (R 2 ), i.e., a probability 
measure in R 2 with finite second moment. Here the ( Wf)t>o are n independent standard Brownian 
motions in R. This result was stated in |5j. In that paper, the authors use a stopping in the n- 
voltage variables which requires finely controlling all trajectories. We prove here a simpler version 
of the result based on a-priori estimates. 


Lemma A.l. Let /o e 


') be a probability with finite second moment, and a set of random 


variables (XLVrt) with law fn. Then (A.l) admits a path wise unique qlobal solution with initial 
datum for l<i<N ' 


Proof. The system (A.l I can be written in 


dZf = / 


2N as the SDE 
dBf+ b(Zf)dt, 
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Figure 1. Permanent (non-transient) regimes of the FhN particle system for 
N = 2 000. Top row: J = 0.1 (A) and J = 3 (B), bottom row: J = 1. The unique 
stationary solution in the small coupling limit analyzed in the manuscript visits 
both attractors transiently (A), while in the high coupling regime (B), the system 
remains around one of the attractors (the system has at least two such solutions). 
In an intermediate regime, the system shows periodic oscillations (bottom row). 


where = (xl,vj,... ,x^ ,v^), a N is a constant 2 TV x 2 N sparse matrix, (B( v ) t > 0 is a standard 
Brownian motion on R 2Ar , and b : U. 2N -*■ U. 2N is a function defined in the obvious way. It is easy 
to see that b is a locally Lipschitz function, moreover, letting (•,■) and |i • || the scalar product and 
the Euclidean norm on R 2Ar respectively, then for all Z N = (a; 1 , *; 1 ,... ,x N ,v N ), 


(Z JV ,b(Z JV )) 


N N t N 

Y l x i (-ax i + bv i ) + Y l v i (v i (v i -X)(l-vi)-x i + I 0 )+j- £ v^v'-v*) 

2 = 1 2 = 1 ^ 2.j=l 


< 


N 

5 > 
2 = 1 


l)a 


N 

2=1 


, 2|2 


J | 2 ^ 


< ^(l + llz^ll 2 ). 


J_ 

N 


N 

Y, v i v j + CN 


i.j=l 


This is a sufficient condition for global existence and pathwise uniqueness (see e.g. j2Sj). □ 


Mean-Field limit. Now we turn to the propagation of chaos property. We already know the 
existence and uniqueness of the particle system , moreover the nonlinear SDE: 


(A.2) 


dv t = ( v t (v t - A)(l -v t ) -x t + l)dt+ J / (v t - v)dft(x,v)dt + dW t , 

dxt = (-ax + bvt)dt 

ft = law(5 t ,u t ), law(x 0 ,u 0 ) = /o- 


is also well-posed for fo e L 1 (M 2 ) n L 1 logL 1 n 
instance, we can sate the 


'), as a consequence of Theorem 2.2 Then, for 


Theorem A.2. Let /o be a Borel probability measure and (Aq, Vq) for 1 < i < N be N independent 
variables with law fo. Let us assume that the solutions to (A.l) and (A.2) with initial data (X q,Vq) 
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and /o are well defined on [0,T] and such that 

(A.3) sup { f (\x\ 2 + \v\ 2 )df t (x,v)\ <+oo, 

[o,t] 1 Jr2 1 

with f t = law(x\,v l t ) (which actually does not depend on i by exchangeability). Then there exists a 
constant C > 0 such that 

(A-4) E[\xl-xl\ 2 + \vi-vl\ 2 ]<^e ct . 

Proof. We start by writing X\ - x\- x\ and Vf = v\-v\. For notational convenience we drop the 
time dependence subindex and take J = 1. Because x\ and x\ are driven by the same Brownian 
motion, we have that 

dV l = (y l (v l - A)(l - v l ) - v l (v l - A)(l ~v l ) - X*) dt+ — ^ (y l t - dt - f {v l - v ) dfifx, v ) dt 
• N j = i J R 2 

dX = (-aX i + bV i )dt, 

We define a(t ) = E[|X*| 2 + |W| 2 ] which is independent of the label i by symmetry and exchange¬ 
ability of the system. It is not hard to see that 

- 2 fE[\Xf]=E[b\XW l \-a\Xf] < b 2 a(t), 

and 

l±E[\V*\ 2 ] = E[y i (u i (u i -A)(l-u i )-u i (u i -A)(l-u i )-A: i )] 

rV iN r n 

+ E [^- E ( v t- v l) dt -V z J R2 (v z -v)dft(x,v)j =:S 1 + S 2 . 

Estimate for Si: Let us first notice that 

v\v { - A)(l - v l ) - €*(€* - A)(l - v l ) = -(H 3 - |t>‘| 3 ) + (1 + A)(|uf - \vf ) - XV i 

= -V\\vf + v* v 1 + |uf) + (1 + A) V l (\v l \ + |v‘|) - XV\ 

therefore 

Si = E[|W| 2 (-|i/| 2 -vW- \vf + (1 + A)(|tA| + |u j; |) - A)] - E[n*], 

and by consequence there is some constant C > 0 such that 

(A.5) Si < Ca(t). 

Estimate for S - 2 : By definition, it holds 

r ... v i N ■ r I 

S 2 = K[v l (vl-v l )- — ^(v J t -J R2 Vdf t (x,v))] 

= E [IH 2 ] - ^ E E (4 - / R2 « dftfx, t,))]. 

Moreover, by symmetry we know that S 2 does not depend on a particular i, therefore we take i = 1 
to get 

S 2 < E[|^| 2 ] + E(E[|^| 2 ]) 1/2 (E[| f (v{ - f R2 vdf t (x,v))\ 2 ]) 112 . 

j-2 

Now, defining Y J = v J t - / R2 v dftfx, v), for j + k, we find that 

E [Y j Y k ] = e[e|Y j ' | (5 1 ,u 1 )]E[y fc | (z 1 ,*; 1 )]], 

but 

E[Y J | (x 1 ,!; 1 )] = E^v( - J ^ vdf t (x,v )j = 0. 
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Hence, fixing j* e {2,..., N} 


E 


[\^{ v3 t ~ f R2 vd M x ’V))\~] = (N- l)E[|nf - ^ vd/ t (x,i;))| 2 ] 

= ( N ~ l ) [ ( w ~ [ vdf t (x,v )) df t fy,w) < C(N - 1), 

Jr 2 ' Jr 2 2 

since the second moment of f t is uniformly bounded in [0, T], Finally we conclude that 


(A.6) 


C 


S 2 < aft) + aft) 1 


Finally, going back to the bounds on a(t), we put together (A.5) and (A.61 to find 


< Caft) + 2a(f) 1/2 ^= < Ca(t) + 


and using Gronwal’s Lemma, 


which finishes the proof. 


t(f) < (a(0) + 


e ct = —e Ct 
N 


□ 


Appendix B. Strong maximum principle for the linearized operator 


In this final appendix we shall extend the result provided in Corollary A.20] to our frame¬ 
work. These local positivity estimates are classical in hypoelliptic equations and they are a nec¬ 
essary condition for Theorem |2.3[ Here, our result is time dependant and by consequence more 
general than it is needed in the applications. 

In the sequel, we shall use the notation 

B r (x 0 ,v o) := {0,u) € R 2 ; \v - v 0 \ < r, \x - x 0 \ < r 3 }, 

and come back to the classical notation Vi,» = D x v and d 2 v = A„. Also, we simplify the problem 
by choosing a = b = 1, but the proof can be easily extended to the general case. 


Theorem B.l. Let f(t,x,v ) be a classical nonnegative solution of 

(B.l) - A v f = Aft, x, v)V v f + B( x, v) V x f + C(t, x, v ) / 

at 

in [0,T) x 0, where f l is an open subset ofM 2 , and A,C : [0,T) x M 2 and bounded continuous 
functions and B(x,v) = x - v. Let (a; 0 , Po) 6 ^ and A and C upper bounds of respectively |A||x,~ 
and |C f |x,~. 

Then, for any r, r > 0 there are constants A, K > 0, only depending on A, C and r 2 lT such 
that the following holds: If B\ r (x o,Po) c H, t < min(l/2, - log(r 3 /2|a;o - Pol)) and f > 6 > 0 in 
[r/2 ,t) x B r (x o,vo), then f > KS in [r/2,r) x B 2r {x o,Po)- 


Theorem |B. 1 1 implies, via covering arguments in variables t,x,v the 


Corollary B.2. If f > 0 solves (B.l I 


m [0,T) x fl and f > 5 > 0 in [0,T) x B r {x o,Po), then for 
any compact set A'c!l containing (xo ,Vq) and for any to e (0, T), we have f > S' > 0 in [to> T) x K 
where S' only depends on A, C, K, H , Xq, vq ,r, to, S. 


Proof of Theorem \B.l\ We only explain how to adapt the proof of Theorem A. 19 given in m\■ Let 
g = e^ t f{t,x, v)-, then g > f and Cg > 0 in (0,T) x LI, where 

L = d t + (v - x) Vx - A„ - Aft, x, v) . 

Next, we construct a particular subsolution for C. In the sequel, B r stands for B r {xo,vo) and 
we define A' t (aio,Po) = vq + (xo - vo)e~ t . 


Step 1. Construction of the subsolution. 
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For t 6 (0, r] and ( x , v) e \ B r let 

P(t,x,v) = a -—- —(v-v 0 ){x-X t ) + 'y -—-, 

with a,/ 3,7 > 0 to be chosen later on. Let further define 

<p(t,x,v) = 6e-» p{t ’ x ’ v) -e, 

where fJ,,£> 0 will also be chosen later on. If we assume that ft 2 < a 7 , then P is a positive quadratic 
form in the variables v - Vq and x - X t . Clearly 

C<p = - f j,5e~ IJ - p £(P), 

where 

£{P) = d t P + (v-x) V X P - & V P + n | S7 V P\ 2 - A(t , x, v ) S7 V P. 

By straightforward computation we find that £ = £\ + £ 2 , with 

UP) = 


t 3 




and 


c( m a (v-v 0 )(x-X t ) 1 

£ A P ) = P - 1 2 - a - t 


( x-X t ) 2 A(t,a;,n)(t;-no) A(t,a;,u)(a;-X t ) 

7 ---a- 1 - + P — 


t 3 t 

Now we notice that £\ is defined by the quadratic form 


t 2 


M q = 


HoP-^-P p+^-nctp 

P+^-fiotP nP 2 -^- 


which is nothing but a quadratic polynomial on (v - vo)/t and (x - X t )/t 2 . As ^ -»■ 00 


ftr M q = fi(a 2 + p 2 ) + 0(1) 


[det M q =l- l [ ' ia f +a/? 7 -/ 3 3 -^-^] + 0(l) 


3 a/3 2 

~2 ^ r 1 2 
both positive quantities if p > a and aj > p 2 . In particular, for f3 = 2a and 7 = 8 a, 

f tr M q =5 a 2 n +0(1) 
jdetM 9 = 2 a 3 ^t + 0(l), 

and letting /i -*• 00 the eigenvalues of M q are of order /i fi 2 and p. So, for any fixed C > 0 we may 
choose a, P, 7 and fi such that 

£ l(P) >_ C^-^l 


Second, if t 6 (0,1) then 


up) >- _ MfcTA. 2/Ji ,. A, 


and making r < 1 , we get, 


£(P) > const y [c( 


P\^(( v ~ v o ) 2 , (x- x ty 


h 


t t 3 

with C arbitrarily large. 

Let us briefly describe the rest of the proof. Recall that (a;, v) $ B r so 
(1) either |v - uq| > r, then £{P) > const. (/3/t)[Cr 2 /t - 1 ], which is positive for C > T/r 2 -, 
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( 2 ) or \x - Xq\ > r 3 , and then, if r < \ min(l, - log( then for any t € [ 0 ,r) 

\X t - Xq\ < r 3 /2 and 


\x-X t \ 2 \x-x 0 \ 2 \X t -x 0 \ 2 r 6 


t 2 2t 2 t 2 4r 2 ’ 

so £ (P) > const. (f3/t)[Cr 6 /At 3 - 1], which is positive as soon as C > Ar 3 /r 6 . 

Summarizing: under the assumptions, we can always choose constants 7 > (3 > a > 1 and 
ay > /3 2 , depending only on A and r 2 /r, so that 

£</?>0, in [0 ,t) x (B Xr ' B r ), 

as soon as r < min(l/ 2 , - log(r 3 / 2 |:ro - uq|)). 


Step 2. Boundary conditions. We now wish to prove that ip < g for t = 0 and for any (x, v ) € 
d(B\ r \ Br)\ then classical maximum principle will do the rest. 

Let us first notice that the boundary condition at t = 0 is obvious {p can be extended by 
continuity by 0 at the initial time). The condition at dB r is also true since V (x,v) e dB r : 
p<S<g. 

It remains to fix the remaining parameters in order to conclude that p < g in dB\ r . From the 
choice of a, /? and 7 , it is easy to see that for any ( x , v ) e dB\ r : 


P(t,x,v) > - 


a/(v-v 0 ) 2 (x-X t y 


t t 3 

notice that we are imposing A > 1. Choosing 


\ a . /, 
J > — mm - 


’ 4r 3 


16 


/r 2 r° \ 


£ = 


5 exp ^ 


pa \ 2 

16 


/r 2 r 6 \\ 
L 1 T ’ r 3 j j’ 


we get ip = S e _ £ < 0 on dB\ r . By consequence p < g on the whole set B\ r . 

Let us finally notice that at this point we have uniform bounds for g on l? 2 r N B r for any 
t e [t/2,t). Indeed, 


P(t,x,v ) < 2 7 1 


(^-^o ) 2 + (x- X t ) 2 


„ / 8 r 2 1026 r 6 \ / r~ r" \ 

< 2 y^-+--— J < 2068 7 max —, — J 


2 . 


t 3 


Then, for A big enough we find I\q > 0 such that 

p(t,x,v) > 5^ exp ^ - 2068 p 7 max ^ —, — j j - exp ^ 


pa \ 2 . / r 2 r 6 \\i 


because 7 = 80 , to find such A it suffices that 

/ r 2 r 3 \ / v 2 r 3 \ 

2068 x 16 x 8 max ^ —, — J < A 2 min —, — J, 


by consequence A depends only on r 2 /r. 

Finally, we find K, A > 0 depending on A, C and r 2 /r such that 


f>K 0 5e TC on [r/ 2 ,r) x (B 2r x B r ). 


□ 


Remark B.3. Let us notice that we can extend Th,eorem \B . 1\ to some cases when A or C are not 
necessarily bounded and £1 = R 2 . It suffices to take any r, r > 0 and fix X (which as we saw only 
depends on a numerical constant and the ratio r 2 It). We can then fix R > 0 big enough, in order 
to have that Xr < R and study the equation into Br, where by continuity A and C attain their 
maximum in the compact set [0,r] x Br. 


















32 


S. MISCHLER, C. QUININAO, J. TOUBOUL 


References 

[1] Abbott, L., and Van Vreeswijk, C. Asynchronous states in networks of pulse-coupled neuron. Phys. Rev 
48 (1993), 1483-1490. 

[2] Amari, S. Characteristics of random nets of analog neuron-like elements. Syst. Man Cybernet. SMC-2 (1972). 

[3] Amari, S.-I. Dynamics of pattern formation in lateral-inhibition type neural fields. Biological Cybernetics 27, 
2 (June 1977), 77-87. 

[4] Amit, D., and Brunel, N. Model of global spontaneous activity and local structured delay activity during 
delay periods in the cerebral cortex. Cerebral Cortex 7 (1997), 237-252. 

[5] Baladron, J., Fasoli, D., Faugeras, O., and Touboul, J. Mean field description of and propagation 
of chaos in recurrent multipopulation networks of hodgkin-huxley and fitzhugh-nagumo neurons. Journal of 
Mathematical Neuroscience 2, 10 (2012). arXiv: 1110.4294. 

[6] Bressloff, P. Stochastic neural field theory and the system-size expansion. SIAM J. on Applied Mathematics 
70 (2009), 1488-1521. 

[7] Bressloff, P. Spatiotemporal dynamics of continuum neural fields. Journal of Physics A: Mathematical and 
Theoretical 45 (2012), 033001. 

[8] Brezis, H. Analyse fonctionnelle. Theorie et applications. Masson, 1983. 

[9] Brunel, N., and Hakim, V. Fast global oscillations in networks of integrate-and-fire neurons with low firing 
rates. Neural Computation 11 (1999), 1621-1671. 

[10] Buice, M., and Cowan, J. Field-theoretic approach to fluctuation effects in neural networks. Physical Review 
E 75, 5 (2007). 

[11] Caceres, M., Carrillo, J. A., and Perthame, B. Analysis of nonlinear noisy integrate and fire neuron 
models: blow-up and steady states. Journal of Mathematical Neuroscience 1 (2011). 

[12] Cai, D., Tao, L., Shelley, M., and McLaughlin, D. An effective kinetic representation of fluctuation-driven 
neuronal networks with application to simple and complex cells in visual cortex. Proceedings of the National 
Academy of Sciences 101, 20 (2004), 7757-7762. 

[13] El Boustani, S., and Destexhe, A. A master equation formalism for macroscopic modeling of asynchronous 
irregular activity states. Neural computation 21, 1 (2009), 46-100. 

[14] Ermentrout, G. B., and Cowan, J. D. A mathematical theory of visual hallucination patterns. Biological 
cybernetics 34, 3 (1979), 137-150. 

[15] Escobedo, M., Mischler, S., and Rodriguez Ricard, M. On self-similarity and stationary problem for 
fragmentation and coagulation models. Ann. Inst. H. Poincare Anal. Non Lineaire 22, 1 (2005), 99-125. 

[16] Ethier, S., and Kurtz, T. Markov processes, characterization and convergence. NY: John Willey and Sons 
9 (1986). 

[17] FitzHugh, R. Mathematical models of threshold phenomena in the nerve membrane. Bulletin of Mathematical 
Biology 17, 4 (1955), 257-278 0092-8240. 

[18] Gamba, I. M., Panferov, V., AND Villani, C. Upper Maxwellian bounds for the spatially homogeneous 
Boltzmann equation. Arch. Ration. Mech. Anal. 194, 1 (2009), 253—282. 

[19] Gualdani, M. P., Mischler, S., and Mouhot, C. Factorization of non-symmetric operators and exponential 
H- Theorem, hal-00495786. 

[20] Herau, F. Short and long time behavior of the Fokker-Planck equation in a confining potential and applica¬ 
tions. J. Fund. Anal. 244-> 1 (2007), 95-118. 

[21] Hodgkin, A., and Huxley, A. A quantitative description of membrane current and its application to con¬ 
duction and excitation in nerve. Journal of Physiology 117 (1952), 500-544. 

[22] Kandel, E., Schwartz, J., and Jessel, T. Principles of Neural Science, 4th ed. McGraw-Hill, 2000. 

[23] Kato, T. Perturbation theory for linear operators. Classics in Mathematics. Springer-Verlag, Berlin, 1995. 
Reprint of the 1980 edition. 

[24] Ly, C., and Tranchina, D. Critical analysis of dimension reduction by a moment closure method in a 
population density approach to neural network modeling. Neural Computation 19, 8 (Aug 2007), 2032-2092. 

[25] Mao, X. Stochastic differential equations and applications. Elsevier, 2007. 

[26] Mischler, S., and Mouhot, C. Exponential stability of slowly decaying solutions to the kinetic fokker-planck 
equation, work in progress. 

[27] Mischler, S., and Mouhot, C. Stability, convergence to self-similarity and elastic limit for the Boltzmann 
equation for inelastic hard spheres. Comm. Math. Phys. 288, 2 (2009), 431-502. 

[28] Mischler, S., and Scher, J. Semigroup spectral analysis and growth-fragmentation equation, hal-00877870. 

[29] Nagumo, J., Arimoto, S., and Yoshizawa, S. An active pulse transmission line simulating nerve axon. 
Proceedings of the IRE 50, 10 (1962), 2061-2070. 

[30] Pakdaman, K., Perthame, B., and Salort, D. Dynamics of a structured neuron population. Nonlinearity 
23 (2010). 

[31] Pakdaman, K., Perthame, B., and Salort, D. Relaxation and self-sustained oscillations in the time elapsed 
neuron network model. SIAM Journal on Applied Mathematics 73, 3 (2013), 1260-1279. 

[32] Pakdaman, K., Perthame, B., Salort, D., et al. Adaptation and fatigue model for neuron networks and 
large time asymptotics in a nonlinear fragmentation equation. 



KINETIC FITZHUGH-NAGUMO SYSTEM 


33 


[33] Touboul, J. On the dynamics of mean-field equations for stochastic neural fields with delays. Physica D: 
Nonlinear Phenomena 2fl, 15 (2012), 1223—1244. 

[34] TOUBOUL, .J., and Ermentrout, G. B. Finite-size and correlation-induced effects in mean-field dynamics. 
Journal of Computational Neuroscience 31, 3 (2011), 453-484. 

[35] Touboul, J., Hermann, G., and Faugeras, O. Noise-induced behaviors in neural mean field dynamics. 
SIAM J. on Dynamical Systems 11, 49—81 (2011). 

[36] Tristani, I. Boltzmann equation for granular media with thermal force in a weakly inhomogeneous setting. 
hal-00906770. 

[37] Villani, C. Hypocoercivity. Mem. Amer. Math. Soc. 202, 950 (2009), iv+141. 

[38] Wilson, H., and Cowan, J. Excitatory and inhibitory interactions in localized populations of model neurons. 
Biophys. J. 12 (1972), 1-24. 

[39] Wilson, H., and Cowan, J. A mathematical theory of the functional dynamics of cortical and thalamic 
nervous tissue. Biological Cybernetics 13, 2 (Sept. 1973), 55—80. 

Stephane Mischler 

Universite Paris-Dauphine & IUF 
CEREMADE, UMR CNRS 7534 

Place du Marechal de Lattre de Tassigny 75775, Paris Cedex 16 

FRANCE 

E-MAIL: mischlerSceremade.dauphine.fr 

Cristobal Quininao 

Universite Pierre et Marie Curie 
Laboratoire Jacques-Louis Lions, CNRS UMR 7598 
4 place de Jussieu F-75005, Paris 
FRANCE 

and Mathematical Neuroscience Team, CIRB 
College de France 

E-MAIL: cristobal.quininao@college-de-franee.fr 

Jonathan Touboul 

College de France Mathematical Neuroscience Team CIRB 
11 place Marcelin-Berthelot 75005, Paris 
FRANCE 

and INRIA Paris-Rocquencourt, Mycenae Team 
E-MAIL: jonathan.touboul@college-de-franee.fr 



